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Abstract 

/ 

A  model  is  developed  for  the  prediction  of  the  seismo-acoustic  noise  spectrum 
in  the  microseism  peak  region  (0.1  to  0.7  Hz).  The  model  uses  a  theory  devel¬ 
oped  by  Cato  [J.  Acoust.  Soc.  Am.,  89 , 1096- 1112 v(  1991  )]Tor  an  infinite  depth 
ocean  in  which  the  surface  orbital  motion  caused  by  gravity  waves  may  produce 
acoustic  waves  at  twice  the  gravity  wave  frequency.  Using  directional  wave  spec¬ 
tra  as  inputs,  acoustic  source  levels  are  computed  and  incorporated  into  a  more 
realistic  environment  consisting  of  a  horizontally  stratified  ocean  with  an  elastic 
bottom.  Noise  predictions  are  made  using  directional  wave  spectra  obtained  from 
the  S  WADE  surface  buoys  moored  off  the  coast  of  Virginia  and  the  SAFARI  sound 
propagation  code,  with  a  bottom  model  derived  using  wave  speeds  measured  in 
the  EDGE  deep  seismic  reflection  survey.  The  predictions  are  analyzed  for  noise 
level  variations  with  frequency,  wave  height,  wind  direction,  and  receiver  depth. 
These  predictions  are  compared  to  noise  measurements  made  in  ECONOMEX  us¬ 
ing  near-bottom  receivers  located  close  to  the  surface  buoys.  Good  agreement  is 
found  between  the  predictions  and  observations  under  a  variety  of  environmental 
conditions. 
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Chapter  1 


Introduction 


Researchers  have  long  been  aware  of  a  peak  in  both  the  seismic  ambient  ground 
motion  spectrum  and  the  ocean  bottom  pressure  spectrum  which  occurs  with  a 
period  of  about  three  to  four  seconds.  This  peak  has  historically  been  called  the 
microseism  peak.  It  has  been  noted  that  the  frequency  of  this  peak  seems  to  occur 
at  twice  the  frequency  of  the  peak  in  the  surface  gravity  wave  spectrum.  Longuet- 
Higgins  [1]  was  the  first  to  develop  a  comprehensive  theory  which  proposed  the 
interaction  of  opposing  surface  waves,  which  closely  approximate  standing  waves, 
as  the  source  of  the  microseism  peak.  He  showed  that  while  the  first-order  pressure 
fluctuations  caused  by  the  gravity  waves  exhibit  an  exponential  decay  in  depth, 
the  second-order  fluctuations  caused  by  the  nonlinear  interaction  of  two  opposing 
waves  does  not  attenuate.  Brekhovskikh  [2]  was  the  first  to  develop  this  idea  into 
a  prediction  of  noise  in  the  ocean.  He  has  been  followed  by  others,  most  recently 
Kibblewhite  and  Wu  [3],  who  all  used  a  perturbation  expansion  as  the  solution 
to  the  wave  equation  with  the  resulting  second  order  solution  as  the  acoustic 
field.  Cato  [4]  has  developed  a  model  which  does  not  rely  on  a  perturbation 
expansion,  and  does  not  require  that  a  standing  wave  approximation  be  made. 
His  theory  predicts  somewhat  higher  noise  levels  for  a  shallow  receiver  than  those 
using  the  perturbation  expansion,  and  a  different  directionality.  In  contrast  to 
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these  predictions,  which  have  the  wave  interaction  as  the  source  mechanism,  are 
those  of  Guo  [5]  who  has  proposed  the  wind  turbulence  acting  on  the  surface 
directly  as  the  source  of  noise  in  the  microseism  peak  region. 

Previous  measurements  of  noise  in  the  deep  ocean  have  been  few,  and  the  re¬ 
sults  have  been  broadly  consistent  with  all  the  above  theories  [6, 7, 8, 9].  What  has 
been  lacking  is  a  direct  accurate  measurement  of  the  source  field  believed  to  be 
causing  the  noise,  either  the  wave  field  or  the  wind  field.  Without  these  quality 
source  data,  broad  assumptions  concerning  the  source  field  must  be  made,  and 
the  potential  variance  in  the  predictions  based  on  these  assumptions  are  typically 
greater  than  the  differences  in  levels  predicted  by  the  different  theories.  Thus 
without  the  accurate  measurement  of  the  source  field,  it  is  difficult  to  evaluate 
the  adequacy  or  inadequacy  of  particular  theories.  In  the  case  of  the  surface 
wave  field,  the  required  measurements  consist  of  directional  surface  wave  spectra. 
This  lack  of  high-quality,  simultaneously  measured,  noise  and  surface  wave  data 
was  the  primary  reason  for  the  deployment  of  the  ECO NOMEX  and  SAMSON 
experiments  [10,11].  An  additional  desirable  feature  of  the  ECONOMEX  experi¬ 
ment  was  its  long-term  nature,  allowing  noise  measurements  under  a  wide  variety 
of  environmental  conditions.  The  measurements  consisted  of  near-seafloor  water 
pressure  and  ground  motion  at  various  depths  from  ~100  meters  to  2500  meters 
in  the  Atlantic  Ocean.  During  the  time  the  ECONOMEX  instruments  were  de¬ 
ployed,  instruments  in  the  same  area  from  the  SWADE  program  were  measuring 
surface  parameters  including  directional  wave  spectra,  wind  speed,  and  other  me¬ 
teorological  quantities.  These  two  coupled  data  sets  provide  researchers  with  the 
measurements  needed  to  test  the  theories  mentioned  above. 

This  work  will  examine  the  model  developed  by  Cato  [12]  for  an  infinite  depth 
ocean  and  will  present  the  derivation  of  that  model  in  detail.  Using  the  SWADE 
directional  wave  spectra  as  inputs,  we  will  estimate  the  resulting  acoustic  source 
level.  We  will  then  incorporate  this  source  level  into  a  more  realistic  environ- 
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ment  consisting  of  a  horizontally  stratified  ocean  with  an  elastic  bottom  using  the 
SAFARI  program  [13,14].  The  geoacoustic  model  for  our  bottom  comes  in  part 
from  the  EDGE  deep  seismic  reflection  study  [15].  The  resulting  ambient  noise 
predictions  will  then  be  compared  to  the  measured  noise  from  ECONOMEX  in 
an  effort  to  judge  the  adequacy  of  our  model. 

In  Chap.  2,  we  review  the  equations  of  motion  and  derive  the  inhomogeneous 
wave  equation  in  terms  of  the  Lighthill  analogy  [16].  The  inhomogeneous  wave 
equation  is  then  solved  in  the  presence  of  a  sound  speed  and  density  discontinuity 
in  Chap.  3.  In  Chap.  4,  we  make  some  simplifications  to  allow  us  to  apply 
the  theory  to  the  motion  of  surface  gravity  waves.  Additionally  in  this  chapter 
we  follow  the  work  of  Cato  [12]  in  deriving  an  expression  for  the  noise  level,  in 
terms  of  the  directional  wave  spectrum,  a  problem  for  which  we  also  introduce 
an  elastic  bottom.  In  Chap.  5  we  describe  the  methods  used  to  make  predictions 
of  noise  levels  using  the  equations  in  Chap.  4,  the  SWADE  directional  spectra, 
and  the  SAFARI  program.  Finally  in  Chap.  7  we  present  the  ECONOMEX  noise 
observations  and  compare  them  to  our  predictions  which,  for  the  most  part,  agree 
very  well.  Our  conclusions  and  recommendations  for  future  research  are  presented 
in  Chap.  8. 
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Chapter  2 


Basic  Equations 


In  any  consideration  of  sound  generation,  the  most  basic  question  to  ask  is  “What 
is  sound?”  That  is,  how  do  we  distinguish  particle  motions  associated  with  acous¬ 
tic  energy  from  that  associated  with  other  types  of  energy  such  as  vortical  energy 
(that  due  to  turbulence)  or  thermal  energy?  In  general,  this  is  a  very  difficult 
question,  for  the  equations  governing  the  three  types  of  motion  are  a  set  of  cou¬ 
pled,  inhomogeneous,  non-linear  partial  differential  equations,  first  derived  in  a 
linear  form  by  Rayleigh  [17],  and  later  in  a  more  general  form  by  Doak  [18].  Doak 
further  showed  that  for  small  Stokes  number  the  set  decouples  into  three  separate 
equations  for  the  acoustical,  vortical,  and  thermal  type  motions.  For  our  model 
of  sound  generation,  we  will  consider  only  this  decoupled,  small  Stokes  number 
approximation,  and  introduce  any  coupling  between  the  three  types  of  motion  via 
external  forcing,  heat  addition,  or  boundary  conditions.  The  acoustic  motion  is 
then  easily  recognized  as  the  small  scale  pressure  perturbation  associated  only 
with  the  internal  elastic  properties  of  the  medium. 

Let  us  consider  the  validity  of  the  small  Stokes  number  approximation  in  the 
case  of  low-frequency  sound  in  the  ocean.  If  we  limit  our  interest  to  a  maximum 
frequency  of  50  Hz,  then  the  maximum  Stokes  number  (u>v/c2)  is  «  10~7,  where  u 
is  the  acoustic  frequency,  u  is  the  medium  viscosity,  and  c  is  the  sound  speed.  We 
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are  then  justified  in  using  the  assumption  of  small  Stokes  number  and  can  take 
advantage  of  the  great  mathematical  and  conceptual  simplifications  it  allows. 

In  our  derivation  of  the  sound  generated  by  a  moving  fluid,  we  will  follow  the 
acoustic  analogy  method  of  Lighthill  [16],  but  will  use  the  pressure  perturbation 
as  our  field  variable  in  preference  to  fluctuations  in  mass  density  as  suggested  by 
Doak  [18].  We  will  first  derive  an  expression  for  the  pressure  perturbations  in  an 
ideal  acoustic  medium  at  rest,  and  then  derive  a  similar  expression  for  a  real  fluid 
which  may  contain  variations  in  density,  variations  in  sound  speed,  and  which  can 
have  existing  within  it  any  type  of  motion.  When  we  compare  the  two  expressions, 
the  difference  will  contain  the  effects  of  motion  in  the  re<d  fluid  which  generate 
sound. 

The  exact  equations  of  mass  conservation  and  momentum  can  be  written  as 


dp  <HpUi)  =Q 

dt  dxi 


(2.1) 


g(guQ  +  a(,  «■■»,■  +  *«)  /i  =  0  (2.2) 

-■t  OXj 

where  p  is  the  mas;,  density,  tq-  is  the  fluid  velocity  in  the  z,-  direction,  cr.y  is  the 
stress  tensor  consisting  of  hydrostatic  and  viscous  stresses,  and  /,•  is  the  body 
force  in  the  Xf  direction.  The  indices  i  and  j  may  take  any  value  1,  2,  or  3, 
corresponding  to  the  Cartesian  coordinate  axes,  and  a  repeated  index  ir.  z  t°rm 
indicates  the  term  is  to  be  summed  over  all  values  of  the  index. 

If  we  consider  the  “acoustic  approximation”,  p  —  p0  =  c0  2(p  -  p0),  where  p  Is 
the  fluid  pressure,  then  Eq.  (2.1)  can  be  rewritten 


1  dp  d(puj) 

cq2  dt  dxi 


(2.3) 


Doak  has  shown  that  this  approximation  is  exact  for  an  inviscid,  non-heat  con¬ 
ducting  fluid,  and  valid  for  a  Stokesian  fluid  to  first  order  in  the  Stokes  number. 
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Equation  (2.2)  can  be  applied  to  our  ideal  acoustic  medium  by  linearizing  it  and 
neglecting  viscous  terms,  leaving 


(2.4) 


(2.5) 


=0. 

dt  dxi 

We  can  now  combine  Eqs.  (2.3)  and  (2.4)  by  eliminating  pm  to  yield 

dt2-  °dx 2 

This,  then,  is  the  familiar  homogeneous  wave  equation  governing  the  pressure 
perturbations  in  an  ideal  acoustic  medium. 

Now  consider  the  exact  equation  of  momentum  balance,  Eq.  (2.2).  This  can 
be  expanded  as 

/it /  -t#  •  4-  A 

(2.6) 


dp  d(pUiUj  +  S(,) 

-ar  +  d^  = - a7, - +  pf" 


where  S  is  the  viscous  stress  tensor,  and  the  hydrostatic  term  has  been  explicitly 
stated.  Next  we  again  eliminate  pm  by  using  the  exact  equation  of  mass  conser¬ 
vation,  Eq.  (2.1).  By  adding  d2p/dx 2  to  both  sides  to  allow  direct  comparison 
with  Eq.  (2.5),  we  arrive  at 


d2p  id2p  _  2d2(pmuj  +Sij)  ,  d2{p-c\p) 
dt 2  c°dx2  c°  dxiXj  dt2 


-  C. 


zdjpfi) 

dxi 


(2.7) 


This  is  the  inhomogeneous  wave  equation  as  derived  by  Doak  [18].  It  is  iden¬ 
tical  to  that  derived  by  Cato  [4]  with  the  exception  of  the  final  term,  which  he 
neglected.  The  pressure  fluctuations  in  a  real  fluid  are  exactly  those  which  would 
occur  in  a  uniform  acoustic  fluid  subject  to  the  external  stress  system  given  by  the 
right  hand  side  of  Eq.  (2.7).  As  c§  is  assumed  constant,  not  only  the  sound  gener¬ 
ation  but  also  the  propagation  effects  of  a  real  fluid  are  included  in  the  equivalent 
stress  system.  It  will  be  our  task  not  only  to  solve  Eq.  (2.7)  in  an  ocean  geome¬ 
try  with  a  realistic  stress  imposed,  but  also  to  separate  out  the  sound  generation 
terms  from  the  propagation  terms. 

It  is  useful  at  this  point  to  consider  the  importance  of  the  term  c\d{pfi)jdxi, 
which  was  neglected  by  Cato  in  his  theory  of  noise  generated  by  surface  orbital 
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motion  [4].  In  the  ocean,  this  term  reduces  to  clg(dp/dz).  We  will  consider 
motions  with  a  periodic  time  dependence  (p  a  e-*"*)  and  a  characteristic  length 
scale  L.  Then  we  can  form  the  ratio  of  the  last  to  the  first  of  the  terms  on  the 
right  hand  side  of  Eq.  (2.7)  as  ( gp/L)/(pu2/L 2).  From  Eq.  (2.1),  we  can  see 
\pU / L\  ~  \ujp\  so  we  can  rewrite  our  ratio  as  g/ Lcj2.  We  can  also  see  from  Eq.  (2.1) 
that  the  first  two  terms  on  the  right  hand  side  of  Eq.  (2.7)  are  of  the  same  order.  If 
we  consider  as  our  worst  case  a  minimum  frequency  of  0.1  Hz,  we  can  neglect  the 
gravity  dependent  term  if  L  25m.  If  we  were  concerned  with  purely  acoustic 
motion,  the  appropriate  length  scale  would  be  the  wavelength  L  =  X  =  c/ f  ~ 
15,000  m.  However  the  length  scale  appropriate  to  the  equivalent  sources  resulting 
from  surface  gravity  waves  would  be  the  wavelength  obtained  from  the  deep  water 
low  frequency  surface  gravity  wave  dispersion  relation  uj 2  =  2-ng/X  which  yields 
L  =  A  ~  160  m.  In  both  cases  we  can  therefore  neglect  the  last  term  in  Eq.  (2.7). 
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Chapter  3 


Solution  to  the  Inhomogeneous 
Wave  Equation 


The  general  solution  to  the  inhomogeneous  wave  equation  is  well  known  (cf.  Strat¬ 
ton  [19]).  For  the  source  term  in  Eq.  (2.7),  it  is  given  by 


PCM)  ~Po  = 


—f 

4irJ 


d2{pUjUj  +  Sjj)  d2(pfcl-p) 

v  [  dyidyj  dt2 


dy 


J_f  \l  dp  J_5r  J_ 
+4xJ  s  r  dn  +  r2  dn?  +  <&r 


dS{ y), 


(3.1) 


1  dr  dp 
c^r  dn  dt 

where  r  —  |x  —  y|,  ft  is  the  outward  normal  to  the  surface  of  integration,  and 
the  integrands  are  evaluated  at  the  retarded  time  r  =  t  —  r/c.  The  volume  of 
integration  above  must  include  all  possible  regions  where  noise  could  be  generated, 
and  the  surface  integral  will  include  the  boundaries  of  these  regions.  We  will  choose 
our  volume  to  be  a  cylinder  of  infinite  height  above  the  sea  surface,  and  extending 


to  the  ocean  bottom.  The  radius  of  this  cylinder  Y  is  assumed  to  be  large  enough 
to  include  any  regions  of  motion  which  could  contribute  to  p(x,t),  and  thus  the 
sides  of  the  cylinder  do  not  contribute  to  the  surface  integral.  We  next  follow 
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Cato  [4]  in  modeling  the  sea  surface  density  discontinuity  as  a  Heaviside  function 
so  that  the  derivatives  in  Eq.  (3.1)  exist.  The  density  is  given  by 


p  =  H(x  -  C)(pw  -  pa)  +  pa 

where  pw  is  the  density  in  water,  pa  is  the  density  in  air,  f  is  the  value  of  x  at  the 
interface,  and  E  (x  —  f)  is  the  Heaviside  function 

H(y)  =  0,  y  >  0, 

H(y)  =  1,  y  <  0.  (3.2) 

We  can  then  expand  terms  containing  spatial  derivatives  of  p  as  follows: 

d2(pUjUj)  _ 
dyidyi 


d 2  Ql 

5  ( Pw  Pa)v-iUj  +  —  ~  (paUiUj) 

dyidyi  dyidyi 

dH  d  d2 

+2-^—  -r—  (Pw  -  PajUiUi  +  (Pw  -  pa)UiUi———H. 


(3.3) 


dy.dyy  dy.% 

We  can  also  combine  Eqs.  (2.1)  and  (2.2)  to  yield 

d2p  _  d2(pu,jUj  +  Qjj) 
dt2  dy^yi 

which  cam  be  used  with  Eq.  (3.3)  to  eliminate  the  derivatives  of  the  Heaviside 
function  in  Eq.  (3.1)  resulting  in 


(3.4) 


4tt(P(x,  t )  -  p0)  = 
dy 


f  —  \H  a"  o  ’  (P*UiUi)  +  (1  ~  H)~  ( paUiUj )  +  ~ ~ 

Jvr  [  dy^  dyidyi  <9y,<9y,J 

r  dylrr  d2  .  .  „  a2  .  .  a2 a2.p, 

/v  r  +  (1  a)ay,3y>-“'Ui)  +  ayiay;  -  aAg1 

1  dp  1  3r  1  dr  dp 


r  1  dp  1  dr  1  dr  dp 
J  s  r  dn  r2  dn^  ^  c^r  dn 


dS(y). 


(3.5) 


Now  we  can  split  the  volume  integral  into  one  integral  over  the  water  volume 
and  one  integral  over  the  air  volume,  noting  that  H  =  1  in  water  and  H  =  0  in 
water,  and  after  recombining  terms  to  form  d2pjdt 2  we  find 


Mp(x,t)  -po)  = 

f  dy  f  d1  .  d2  .  p  . 

Jv.T  [  W','UiU’  +  +  w{4  ~  pm\ 

r  dy  f  d2  . 

+Jv.T  + 

r  Tl  dp  1  dr  1  dr  dp 
+J  s  [r  dn  **"  r2  dv?  e^r  dn  dt 


)  +  dt^4  p*\ 


dS(  y). 


(3-6) 


As  before,  the  integrands  must  be  evaluated  at  time  r  =  t  —  rfc 

The  next  step  is  to  apply  the  divergence  theorem  in  a  manner  first  described 
by  Curie  [20],  and  presented  in  detail  by  Cato  [12].  This  is  not  a  straightforward 
process,  as  the  terms  to  be  evaluated  depend  on  y  both  directly  and  also  through 
the  dependence  on  r.  Curie  stated  and  Cato  demonstrated  that 


f  a'Fu  dy 

f  dF^  dy 

r  jl 

\dFij  1] 

J  v  dy{dyj  r 

dxjJ 

v  dyj  r 

J  vdyi 

.  dVi  r. 

dy 


=  r  i  dS(y) 
J  s  *  dyj  r 


(3.7) 


where  are  the  direction  cosines  of  the  outward  normal  of  the  surface  enclosing 
V .  We  can  then  repeat  this  operation  a  second  time  to  yield 

dFij  dy  d  f  dy  r  dS{ y) 


f  ^  *f  Fij~  =  / 

J  v  dyj  r  dXiJ  v  r  Js  r 

We  can  now  combine  Eqs.  (3.8)  and  (3.9)  to  achieve  Curie’s  result 

/, 


(3.8) 


d2F{j  dy  _ 
vdy{yj  r 


dy  +  dr  dy  +  J  ^ 

oxidxjj  v  r  dXiJ  s  r  J  s  dyj  r 


dy 


dFij  dS(y) 


(3.9) 
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We  can  also  easily  transform  the  surface  integral  in  Eq.  (3.1)  as  follows.  First  we 
resolve  the  normal  derivatives  into  Cartesian  coordinates  which  yields 


f  \l  dp  1  dr  1  dr  dp] 

J  s[r  dn  + r2  dnP  +  clr  dn  dt  T  ^ 

f  .  \l  dp  1  dr  1  dr  3pl  .  . 

J  si[rlfyi+^l&iP+&d^idi\T 

By  substituting  dr/dy ,•  =  —dr/dxi  and  dp/dy ,•  =  d(pSij)/yj,  we  have 

r  [l  dp  1  dr  1  3r  dp 1  .  . 

/  U  — z—  p  H  2~a  *57  ds(y)  ~ 

J  s  [r  dy{  r2  dy{  cgr  dy{  dt\r 

f  i  a(p6ij)  dS(y)  f  i  f1  dr  1  dr  jg|  ds(y) 

J  s  *  dyj  r  Js'[r  dxP  c§  dxi  dt  r  r 


Next  we  note  that 


so  that  we  can  write 


f  j  d(PSii)  dS(y)  f  ,  f1  d*  ,  1  9r  dp]  dS( y) 
J  s  ‘  dyj  r  J  s  *  [r  dnP  dxi  dt  r  r 

=  / 

J  s  dxi  r 


Combining  Eqs.  (3.10),  (3.11),  and  (3.12)  yields 


f  \l  dp  1  dr  1  dr  dp]  .  . 

Js  [71fa  +  7ilhlr  +  cf'ifa!)i\'  *y'~ 

f  +  f  <**"! M. 

J  s  dyj  r  Js  dxi  r 


(3.10) 


(3.11) 


(3.12) 


(3.13) 


Now  we  substitute  Eqs.  (3.9)  and  (3.13)  into  Eq.  (3.5).  After  noting  that  the 
source  motions  at  infinity  do  not  contribute,  we  can  separate  the  surface  integrals 
enclosing  the  water  volume  into  one  at  the  air-sea  interface  designated  5,-  and  one 
at  the  ocean  bottom  designated  Si,  and  can  then  combine  the  surface  integral  at 
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the  bottom  derived  by  applying  the  divergence  theorem  to  V.  with  the  surface 
integral  of  Eq.  (3.13).  This  yields 


4jr(p(x,t)  ~Po)  = 

a2 


aXidxJ +  /v.^i(|  ”  "-It 

^/v>“'“'  +  +  /v.|i (|  -  -Jr 


dxidxj 

"/,/*  +  •%)]  ^ 

+d7if sMpaUiUi  +  Si^~~  +  /5/«  +  *%)] 

+dxifsJi(pwUiUi  +  Sii  +  p6h)~~ 

+lstli  [%^‘oU,Uy  +  SV  +  rff)]  (3.14) 

We  can  transform  the  y  derivatives  to  time  derivatives  by  noting  from  Eq. 
(2.2)  that 

+  S»  +  P4i)j  = 

and  by  defining  Ap  =  (p.  _  p.),  we  can  finaUy  write 
4yr(p(x,<)  -po)  = 

+ s*>f + /v.^(|  -  "-It 
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(3.15) 


We  can  now  begin  to  interpret  the  meaning  of  Eq.  (3.15)  in  terms  of  the  sources 
of  classical  acoustics.  The  double  spatial  derivatives  of  the  volume  integrals  are 
recognized  as  distributions  of  quadrupole  sources  in  the  volume  V ,  and  the  single 
spatial  derivatives  of  the  volume  integrals  are  recognized  as  dipole  source  distri¬ 
butions  in  the  volume  V.  The  spatial  derivatives  of  the  surface  integrals  at  the 
interface  are  seen  as  dipole  sources  caused  by  the  motion  of  the  surface  of  the 
density  discontinuity,  and  similarly  the  surface  integrals  of  the  time  derivatives 
are  seen  as  monopole  source  distributions  caused  by  the  motion  of  the  surface  of 
the  density  discontinuity.  These  interpretations  are  those  given  by  Lighthill  [16]. 

In  order  to  interpret  the  terms  involving  the  double  time  derivatives,  we  must 
make  a  further  approximation  to  allow  us  to  simplify  Eq.  (3.15),  although  it  is 
one  that  is  frequently  made  in  classical  acoustics.  If  we  assume  that  the  fluids  are 
isentropic,  then  the  density  and  pressure  are  simply  related  by  dp/ dp  =  c2,  where 
c  is  the  local  speed  of  sound  [21].  Also,  if  entropy  is  conserved,  the  processes 
are  reversible  and  the  dissipative  effects  of  viscosity  can  be  ignored.  (Note  we 
must  now  treat  the  effects  of  attenuation  due  to  absorption  separately.)  With  no 
viscosity,  <S,/  =  0,  or  ertJ  =  p<5,;- .  We  can  thus  write 


& 

t2 


or  using  Eq.  (3.4)  we  have 


(3.16) 
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We  can  again  apply  the  divergence  theorem  twice  to  Eq.  (3.16)  to  yield 

J  vt2  PJ  r 

/’  /  (1  i 

axidxjJ  v  \cg 

\  dy 

1  (pxiiUj  +  pSij)— 

+dxj s’  (eg  ~  1 

\  dy 

j  (pUjU,  +  pSij)— 

[ldUC2  * 

1  „„  1 

Js'dt[{c2  , 

)  PU«J  r  ’ 

-y/s[u-iH  r  (317) 

From  the  above,  we  note  that  we  can  again  identify  both  a  quadrupole  and  a 
dipole  volume  source  distribution,  and  both  a  dipole  and  a  monopole  surface 
distribution.  Next  we  combine  Eqs.  (3.17)  and  (3.15)  resulting  in 

47r(p(x,<)  -p0)  = 

d2  f  [c2,  .  .  dy 

v.  +  p6ii)  ~  P^J  T 

sli  ^pwc2'°  ~  PoC^u,Uj + _  c*)p<5‘>] 

s^ai  [(p^  “  PoC*H  ~r^ 
d  f  t  chr  \  c  \ds[y) 

St  Cg  r 

This  result  shows  that  the  sound  received  at  a  point  is  the  sum  of  apparent 
sources  distributed  throughout  the  source  volume,  at  the  air-sea  interface,  and  at 
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the  ocean  bottom.  The  effects  of  the  real  ocean,  such  as  refraction,  are  included 
in  Eq.  (3.18).  The  reflection  from  the  ocean  surface  and  bottom  are  also  modeled 
in  Eq.  (3.18).  Our  next  step  will  be  to  further  simplify  Eq.  (3.18)  to  identify  the 
real  sources  which  transform  other  types  of  energy  into  acoustic  energy,  and  to 
show  how  the  sound  generated  by  these  real  sources  is  modified  by  a  homogeneous 
ocean. 
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Chapter  4 


Application  to  Orbital  Motion  in 
the  Deep  Ocean 


Equation  (3.18)  is  quite  general  and  could  in  theory  be  used  to  predict  the  acoustic 
pressure  for  a  given  source  mechanism,  but  would  require  a  knowledge  of  the 
particle  motions  u(-  throughout  the  volume  of  the  ocean  and  on  the  boundaries. 
In  addition,  the  effects  of  a  real  fluid  such  as  refraction  are  modeled  as  apparent 
sources,  which  further  confuses  the  issue.  We  will  now  begin  to  make  assumptions 
which  will  greatly  simplify  Eq.  (3.18)  to  allow  us  to  use  it  to  predict  infrasonic 
noise  generated  by  the  sea  surface  orbital  motion. 

We  first  assume  a  plane  wave  sea  surface  elevation 


f  =  acos(#e  "X1  —  at), 


(4.1) 


where  f  is  the  sea  surface  elevation,  a  is  the  amplitude,  k  is  the  wavenumber 
vector,  x*  =  Xjt-f  x 2j  is  the  horizontal  position  vector,  and  a  is  the  radian  gravity 
wave  frequency.  The  solution  to  the  linear  equations  of  motion  yields  a  velocity 
potential  <f>  at  depth  x$  =  z  given  by 


i  <7<xcosh/c(z  +  d)  .  ,  , 

4>  =  -  - .  T — J - sin  K  •  x  -  at), 

k  smh  Kd 


(4.2) 
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where  d  is  the  total  water  depth  and  k  =  |#c|  [22].  Ignoring  surface  tension,  this 
leads  to  the  dispersion  relationship 

a2  =  gKta.nh.Kd,  (4.3) 

where  g  is  the  gravitational  constant.  If  we  restrict  ourselves  to  water  deep  enough 
to  satisfy  icd  >  tt/2,  Eq.  (4.3)  reduces  to  the  deep  water  surface  gravity  wave 
dispersion  relation 

a2  =  Kg.  (4.4) 

We  then  find  the  particle  motions  to  be 

x\  =  ae“°  sin(*c  •  xj,  —  at),  t  =  1,2, 
z  =  ae**°  cos(#c  •  xj,  —  at),  (4.5) 

where  the  initial  position  of  the  particle  is  (Xq,  Zo).  We  can  see  that  the  particles 
describe  a  circular  path,  and  hence  the  term  orbital  motion.  Now  we  must  de¬ 
termine  the  minimum  water  depth  for  which  we  may  neglect  the  orbital  motions 
at  the  bottom.  We  will  also  need  to  find  the  minimum  water  depth  for  which 
Eq.  (4.4)  holds.  If  we  consider  a  minimum  gravity  wave  fi  equency  of  0.05  Hz, 
the  deep  water  dispersion  relationship  is  valid  in  water  depths  >  200m.  To  find 
the  depth  at  which  we  can  negiect  orbital  motion  on  the  bottom,  we  consider  the 
c~K*  decay  with  depth  and  we  find  the  orbital  motions  are  5%  of  the  surface  value 
at  a  water  depth  >  300m.  We  will  therefore  consider  to  be  the  effective  source 
generation  region  to  be  <300  m.  Outside  this  region  we  may  neglect  the  volume 
integral  of  Eq.  (3.18)  and  also  neglect  the  orbital  motions  on  the  bottom. 

We  can  eliminate  any  apparent  sources  caused  by  refraction  by  considering 
the  case  of  an  isovelocity  ocean,  or  cw  =  Co,  and  by  considering  the  air  also  to  be 
isovelocity,  c„  =  t£.  We  can  also  eliminate  a  number  of  terms  in  Eq.  (3.18)  by 
noting  that  c*  »  c*  and  pwc*  »  pac\.  Thus,  we  neglect  the  air  volume  integral, 
and  we  assume  the  upper  half-space  to  be  a  vacuum.  We  have  now  reduced  the 
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motions  on  the  bottom  to  those  due  only  to  the  motion  of  the  acoustic  wave. 
Thus,  in  Eq.  (3.18)  we  can  eliminate  from  the  bottom  integral  the  p«,u,-uy  and 
d(pu) /dt  terms,  since  for  an  acoustic  wave  pu  <C  p.  After  all  these  simplifications, 
we  are  left  with 


t  .  .  .  .  d2  f  dy 

Mp(x,i)-po)  = 

d  f  ,  r  c  ,  dS( y) 

~  Ji  [PvWi  +  Phil 


+  I s!jdt [Pu,Ui^  r  +  dxj s!jP6ij  /  •  (4'6) 
One  can  now  identify  each  term  in  Eq.  (4.6)  and  relate  it  to  a  specific  physical 
mechanism.  The  quadrupole  volume  integral  involving  pu,-uy  is  identified  as  the 
volume  contribution  to  the  second  order  pressure  fluctuation  first  proposed  as 
the  cause  of  microseisms  by  Longuet-Higgins  [l]  .  The  dipole  surface  integral 
involving  pu,uy  is  identified  as  the  sea  surface  contribution  to  this  mechanism. 
The  monopole  surface  integral  containing  ^pu,-  is  identified  as  the  first  order 
pressure  effect  which  attenuates  with  depth.  The  apparent  sources  due  to  the 
surface  and  bottom  reflections  are  contained  in  the  two  p%  terms  .  We  can  see 
in  this  case  that  the  bottom  does  not  influence  the  real  source  terms  directly,  and 
the  real  sources  caused  by  the  transformation  of  mechanical  energy  to  acoustic 
energy  can  be  evaluated  first,  and  the  bottom  effects  added  later. 

Cato  has  shown  [4]  the  quadrupole  contribution  to  be  <  2a /ck  times  the  dipole 
contribution  in  the  far  field,  which  is  defined  by  tk  >  1.  For  a  minimum  wave 
frequency  of  0.05  Hz,  this  is  then  valid  at  depths  greater  than  100  m  below  the 
region  of  effective  volume  generation,  which  we  found  to  be  no  greater  than  300 
m.  Thus  we  are  able  to  ignore  the  quadrupole  term  for  receivers  greater  than  400 
m.  Since  we  are  mainly  interested  in  the  microseism  peak  region,  we  will  also 
neglect  the  monopole  term,  although  it  will  be  significant  at  moderate  depths  for 
frequencies  below  the  microseism  peak. 
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In  finding  the  acoustic  source  level,  we  must  make  two  major  calculations. 
The  first  is  to  relate  the  source  acoustic  power  density  P(u>,z)  to  the  frequency- 
wavenumber  spectrum  of  the  source  motions  pu,uy.  We  must  then  relate  this 
spectrum,  say  ${u>,k),  to  fl(cr)  which  is  the  readily  measured  power  spectrum  of 
the  sea  surface  elevation  f.  Both  of  these  derivations  have  been  performed  by 
Cato  [4,12],  and  we  will  present  the  important  steps  for  the  dipole  sources  below. 
We  must  then  incorporate  these  source  levels  into  an  appropriate  ocean  model  to 
find  the  noise  power  for  an  ocean  receiver. 

We  start  by  redefining  some  terms.  Let 

pD{x,t)  =  ~  I  sljG(r)Wi3(y,T)dy,  (4.7) 

where  W,-y  =  piqu,-,  G(r)  =  l/(47rr),  and  y  is  the  source  coordinate  on  the  sea 
surface.  Since  we  will  perform  the  integration  over  the  mean  sea  surface,  we  can 
immediately  align  the  x3  axis  with  the  depth  axis  z  and  then  li  =  h  =  0,  /3  =  1, 
and  Wij  =  Wi  =  pu,u3.  We  assume  W,  to  be  temporally  stationary  and  ergodic, 
and  spatially  homogeneous.  Then  the  autocorrelation  function  of  pD(x,t)  is 

^p(M  +  ?)  =  ((PD(x,f)p£>(x,t  +  q)),  (4.8) 

and  since  the  process  is  ergodic 

Rp{t  +  <;)  =Rp{<;)  =  rlim  ^  J  ^  P£»(x,f)pn(x,f  +  q)dt,  (4.9) 

and  the  power  spectrum  is 

Qd(u)  =  Um  f  ~f  pD{x,t)pD(x,t  +  q)dteiu,id<; 

T -»oo  J-oo  11  J-T 

=  I'slCwJ-t  [ikfsG{rm*’r)d9. 

X  [^7  Js  GWWtf,  r')#'j  dt  dq.  (4.10) 
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Since  y  and  y'  are  independent,  we  can  define  a  separation  parameter  tj  in  a  cross 
correlation  function  of  W ,• 


Rw„  =lim±[T  W.tfJWW  +  n,t  +  ()dl,  (4.11) 

7— »oo  2,±  J -T 

and  we  can  further  define  the  cross  spectrum  of  Wt- 

Rwi((ri , c)e*u,? ds.  (4.12) 

J—oo 

We  also  note  that 

.  fOO  A\r 

h“)  =  J  ^  (4-13) 

Finally  we  can  combine  Eqs.  (4.10)-(4.13)  and  after  making  some  substitu¬ 
tions,  find  that 

/OO  A\r 

$tI(w,  k)fr,(w,k,  z)H;{w,  k,  (4.14) 

where 

H.KM)  =  («-k*//'<(“r/'+k',T)  •  <4-I5> 

This  completes  the  first  step  in  determining  the  source  levels  from  surface 
orbital  motion.  We  have  the  source  level  QD  as  an  integral  of  the  wavenumber- 
frequency  spectrum  of  the  orbital  motions  times  the  coupling  factors  JJ.IZ/.  Cato 
has  calculated  Hi  analytically  for  an  infinite  ocean  radius,  and  the  values  can  be 
calculated  numerically  for  a  finite  ocean  radius  [4].  Cato  provides  a  thorough 
discussion  of  the  behavior  of  the  coupling  factors  [4j  and  has  shown  the  spatial 
dependence  of  the  coupling  factors  to  be 

Hi  a  exp  j— t  |k  •  x  +  zyj(u2/c2)  —  &2J  j  .  (4-16) 

One  can  see  that  the  spatial  dependence  is  entirely  oscillatory  in  the  horizontal 
direction,  but  is  oscillatory  in  the  vertical  only  for  k  <  u)jc.  For  k  >  uj/c  the 
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coupling  factor  magnitude  has  an  exponential  decay  with  depth.  This  will  lead  to 
an  important  approximation  later  in  the  development. 

Cato  has  also  provided  a  convenient  form  for  the  coupling  factors  in  terms  of 
dimensionless  quantities  M  =  u/kc,  X  =  rfw/c,  X0  =  Rlj/c  and  A  =  zw/c: 

Hi(<jj,k,z)  =  —icoaaH0(<jJtk,z)/M, 

Hi(u,k,z)  =  —ismaHo(u>,k,z)/M, 

Hs{u,k,z)  =  H0{u,k,z ),  (4.17) 

where  a  is  the  angle  between  k  and  the  y\  axis,  and  Hq  is  the  monopole  coupling 
factor  _____ 

1  rx 0  Xe~wx^+A'' 

ffo(w,  k,  z)  =  -  fo  ^  MX/M)  dX,  (4.18) 

J0  being  the  Bessel  function  of  the  first  kind,  order  zero. 

We  must  now  relate  the  frequency-wavenumber  spectrum  4>,j  to  the  surface 
wave  height  power  spectrum  fl.  We  start  by  relating  the  spectrum  of  u,-  to  that 
of  pU{Uj.  Cato  has  shown  [12]  that 

$«yim(w,k)/p2  =  tf,i(w,k)  *  ^,m(w,k)  -f  ¥im(w,k)  *  ¥fi(w,k),  (4.19) 


where  *  denotes  convolution  and  '&,-j(w,k)  is  the  power  spectrum  of  u,-.  We  can 
expand  one  of  these  terms  in  polar  coordinates,  with  a  being  the  gravity  wave 
frequency,  as 


[$,yjm(w,  k,  a)]i  =  p2  f  f  f  ¥,-((<7,/c,7) 

Jk=Q  J 7=0  J <7— — -oo 


7=0  Ja=—oo 


(4.20) 


where  k'  =  k  —  n,  with  7,  7'  being  the  angles  k,  k!  make  with  the  y\  axis.  We  can 
define  a1  =  w  —  <7,  and  can  now  identify  o  and  </  as  the  frequencies  and  k  and  k' 
as  the  wavenumber  vectors  of  the  two  interacting  waves  u,uy,  and  uj  and  k  as  the 
resulting  acoustic  frequency  and  wavenumber.  We  can  relate  the  wavenumbers 


26 


and  angles  by 

k  cos  a  =  k  cos  7  +  k!  cos  7', 

A;  sin  a  =  /csin7  +  tc's'ini'.  (4*21) 

Since  cr  and  k  (and  a1  and  k')  are  associated  with  the  orbital  motion  of  surface 
gravity  waves,  they  are  uniquely  related  by  Eq.  (4.4).  This  is  not  true  for  u>  and 
k.  We  can,  however,  use  the  behavior  of  the  coupling  factors  previously  described 
to  find  approximations  for  u  and  k  in  terms  of  a  and  k. 

We  will  again  follow  Cato  [12]  in  noting  the  exponential  decay  of  the  coupling 
factors  as  k  exceeds  w/c,  and  we  can  choose  a  value  P  such  that  we  will  only 
consider  a  Fourier  component  to  ensonify  a  region  if  k  <  /?w/c.  From  Cato’s 
examination  of  the  coupling  factors  [4],  we  can  choose  /?  =  5  for  a  practical 
receiver  depth.  We  can  now  write 

—  /?w/c  <  |k|  <  /?w/c,  (4.22) 

but  |k|  >  |ic|  —  |jc'|  so  we  may  use  the  deep  water  dispersion  relationship  to  write 

—  Pw/c  <  a2/g  —  (w  —  o)2 /g  <  pu/c,  (4-23) 

or 

jw  —  2o|  <  @g/c.  (4-24) 

Thus,  the  acoustic  frequency  will  be  within  3%  of  twice  the  gravity  wave  frequency. 
We  can  say  with  reasonable  accuracy  that 

u  ~  2a  ~  2 a'.  (4-25) 

We  now  have  formally  obtained  the  frequency  doubling  effect,  where  two  gravity 
waves  of  the  same  frequency  interact  to  generate  an  acoustic  wave  of  twice  that 
frequency.  We  can  similarly  show  that  k  ~  k'. 

Next  we  define  the  two-dimensional  wavenumber  spectrum  of  u,- 

%«{*,!)=  f  *(<^*,7)!p-  (4-26) 

J  —00  Z7T 
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Since  a  and  k  are  related,  we  can  define  a  as  the  specific  value  of  a  determined 
by  Eq.  (4.4).  Also,  since  k  and  ui  are  not  variables  in  the  convolution  integral  Eq. 
(4.20),  /c'  =  |k  —  #c|,  and  we  have  the  two  relations  of  Eq.  (4.21),  we  can  find  k 
and  k'  in  terms  of  7.  We  will  define  the  values  of  0  and  k  thus  determined  as  dn 
and  k .y.  Now  we  will  make  the  assumption  that  the  two-dimensional  wavenumber 
spectrum  is  separable  into  a  one-dimensional  spectrum  times  an  angular  spectrum, 

X./f'c.'y)  =  27rx,i(«)G(a,7)/#c,  (4.27) 


where  the  angular  dependence  satisfies 

/2*G(<7,'y)<f7  =  l.  (4.28) 

Jo 

After  substituting  Eqs.(4.26)  and  (4.27)  into  Eq.  (4.20),  making  the  variable 
substitutions  mentioned  above,  and  also  substituting  d/c  =  (d/c /do)  do,  we  note 
that  we  can  eliminate  two  of  the  three  integrals  in  the  convolution.  Then 

[$.y/m(w,*,a)]i  =  2?rp2  f  Xi»(«)Xim(K,)G(<r,7)G(< 

J'l—O 

1  3k~  , 

x  «;(y)  do^  1 

We  arrive  at  a  similar  expression  for  the  second  term  on  the  left  hand  side  of  Eq. 
(4.19)  by  substituting  subscripts.  We  can  write  the  sound  pressure  spectrum  as 

QdH  =  f  2? rp2  f  (x.i(K)X/m(/c') 

J-  00  J~i=0 


+  X.m(K)xJj('c')]G(<T>7)G!(^',7') 


V  TT*  — 

do  k'  lm( 2tt)2' 


(4.30) 


Cato  has  shown  that  for  significant  ensonification  of  the  noise  field  and 

o  can  vary  by  only  a  small  proportion  of  their  values,  and  thus  can  be  regarded 
as  constants  in  the  integration.  Also,  G(< 7,7)  varies  relatively  slowly,  and  k  ^  k', 
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so  we  can  substitute  G(o,~f)G{o,'i  +  ir)  for  G(a,'i)G(al,'f).  We  can  differentiate 
Eq.  (4.5)  to  find 

Xu  =  cos2  7  •  Xss, 

Xm  =  sin2  7  •  Xs3> 

Xis  =  cos-7  •  xss, 

Xm  =  sin^y-xsa-  (4.31) 

We  further  define  the  one-sided  frequency  spectrum  of  us  as  ft  (a)  and  can  write 

ftW  =  (*),  (4.32) 

and  since  u3  =  df/dt,  we  have  the  frequency  spectrum  of  the  wave  height  in  terms 
of  the  wavenumber  spectrum 

flM  =  ~X»  («)•  (4-33) 

We  can  now  see  that 

Xi»(*)  =  (4.34) 

where  gu  can  be  found  from  Eq.  (4.31)  and  gn  =  1.  We  note  that 


gu(l)9ss{l  +  tt  )$«  =  gai^gisil  +  rfhi, 


(4.35) 


and  so  both  terms  of  Eq.  (4.30)  are  identical.  We  can  then  write  the  sound 
pressure  in  terms  of  fi(<r), 

,64  do-  - -  r*®  r2r _ dadk 


QdH  =  2P2~{l3(a)Iau(&)  JQ  fQ  X  (4.36) 


where 


gu{l)giz{l  +  *)G{o,'i)G{a,7  +  n)di.  (4.37) 

Jo 

We  can  see  that  Iail  +  =  Ias3. 
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Next  we  separate  out  the  a  dependence  of  the  Hi,  perform  the  a  integration 
of  Eq.  (4.36),  and  note  that  Hi  =  H2  to  find 

QdH  =  p2^|^n2(d)/asi(d)  JQk°  (H1H1  +  2  HsHs)  k  dk,  (4.38) 

A  A 

where  Hi  =  Hi[w,k,z).  Finally,  we  use  the  dispersion  relationship  to  substitute 
k!  =  a^/g,  8&/k  =  g/(2o),  and  we  use  Eq.  (4.25)  to  substitute  o  =  o'  =  o//2, 
and  we  have  the  result 

QdH  =  ^-wn2(w/2)/ass(o;/2)  JQ  (HlHi  +  2HsH;)kdk.  (4.39) 

This  gives  us  the  received  pressure  for  an  infinite  depth  ocean  in  terms  of  known  or 
measurable  quantities.  We  can  measure  fi(o')  directly  or  use  an  empirical  spectrum 
based  on  wind  speed.  We  can  similarly  measure  or  imply  from  the  wind  speed 
<2(0,7)  and  thus  calculate  Ia3S  (0) .  HtHf  can  be  calculated  for  a  given  geometry, 
and  p  and  g  are  physical  constants.  We  thus  have  an  equation  which  will  allow  us 
to  predict  the  noise  spectral  levels  for  a  receiver  in  an  infinitely  deep  ocean. 

We  must  now  turn  to  the  problem  of  finding  the  noise  levels  in  a  more  realistic 
ocean.  We  will  consider  the  case  of  horizontally  stratified  media  in  the  water 
column  and  in  the  bottom.  Schmidt  and  Kuperman  have  shown  that  we  can 
write  the  noise  intensity  generated  by  a  horizontal  distribution  of  homogeneous 
sources  at  depth  d  as  [23] 


P(v,z,z')  =  g2(w,z')l£(w,z,z'), 


(4.40) 


where  g2(u,  z1)  is  the  acoustic  monopole  source  strength  of  the  distributed  sources 
of  order  m,  and  we  can  call  the  bottom  gain: 

RTrm  f  00 

Tl(u>,z,z’)  =  J0  |<7(M>*')![ K2{z')  -  k2]k  dk.  (4.41) 

In  this  equation  K(z )  =  w/c(z),  and  g  is  the  depth-dependent  Green’s  function 
satisfying 

g  +  |j = 


(4.42) 


We  can  eliminate  the  depth  dependence  in  z1)  by  normalizing  the  source 
strength  to  that  produced  by  the  same  source  distribution  in  an  infinitely  deep 
ocean,  Q(w),  where  to  first  order  =  Q2(uj)/16n(zt)2.  Now,  provided  z1  is 

chosen  small  enough  compared  to  the  wavelength,  Eq.  (4.40)  is  approximately  in¬ 
dependent  of  z*,  and  we  have  only  to  find  the  equivalent  monopole  source  strength 
and  evaluate  T(u,z,  z1)  to  find  P(u>,  z).  While  the  solution  to  Eq.  (4.41)  for  an 
elastic  bottom  is  not  available  analytically,  there  are  numerical  solutions  available. 
We  can  easily  relate  Qn(w)  to  its  equivalent  monopole  source  strength  by  noting 
from  Eq.  (4.17)  that 

h,h[  +  h,h;  = 

HSH;  =  -i-fliflj.  (4.43) 

c  az 

Thus  the  combined  horizontal  components  are  of  vertical  order  m  =  1  and  the 
vertical  component  is  of  order  m  =  2.  We  can  now  write 

PdH  =  [^-wnI(u/2)J«,,(w/2)j  x  1*°  HoH-JhPkdk 

+  T2(uj,z)  fk°-HoH*kdk  ,  (4.44) 

JO  c 

where  the  index  on  T  indicates  the  order  of  the  vertical  model  used  in  Eq.  (4.41). 

We  now  have  the  full  solution  to  the  received  pressure  spectrum  in  a  horizon¬ 
tally  stratified  ocean.  We  will  next  look  at  the  specific  techniques  and  data  used 
to  generate  predictions  of  the  noise  pressure  spectrum. 
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Chapter  5 


Prediction  Techniques 


5.1  SWADE  data 


The  SWADE  project  (Surface  Wave  Dynamics  Experiment)  was  an  effort  to  char¬ 
acterize  the  sea  surface  using  a  variety  of  sensors  and,  at  the  same  time,  to  measure 
other  relevant  environmental  parameters  [24].  The  project  included  several  pitch 
and  roll  surface  buoys,  satellite  radar  backscatter  measurements,  SWATH  ship  ar¬ 
ray  deployments,  and  aircraft  overflights.  The  long-term  deployment  of  the  pitch 
and  roll  buoys  is  the  element  of  the  experiment  of  direct  use  to  us  in  our  effort  to 
predict  noise  generated  by  surface  orbital  motion. 

SWADE  and  ECONOMEX  instrument  locations  are  shown  in  Fig.  1,  and 
the  overall  experimental  schematic  is  shown  in  Fig.  2.  Deployment  locations  of 
SWADE  instruments  are  listed  in  Table  5.1. 

The  buoys  provided  one  complete  set  of  measurements  including  a  spectral 


|  Instrument  name 

Type 

Latitude 

Longitude 

Water  depth 

DISCUS  E 

Pitch  and  roll 

73°  23. 5  W 

2670m 

DISCUS  C 

Pitch  and  roll 

74°  23. 5 W 

102m 

DISCUS  N 

Pitch  and  roll 

73°  38. 9 W 

115m 

Table  5.1:  Relevant  SWADE  instrument  summary. 
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estimate  each  hour,  allowing  us  to  update  our  estimate  of  the  noise  each  hour. 
The  spectral  data  from  the  SWADE  buoys  was  provided  for  frequencies  from  0.03 
Hz  to  0.34  Hz  in  0.01  Hz  bins,  allowing  acoustic  predictions  from  0.06  Hz  to 
0.68  Hz.  The  SWADE  buoys  which  were  closest  to  the  ECONOMEX  instruments 
were  the  Discus  E  and  Discus  C  buoys.  For  predicting  the  noise  at  the  deep 
ECONOMEX  site,  the  data  of  Discus  E  was  used,  while  predictions  of  noise  at 
the  three  shallow  ECONOMEX  instruments  were  based  on  the  data  of  Discus  C. 
The  SWADE  project  had  deployed  a  SPAR  buoy  near  the  shallow  ECONOMEX 
site  with  better  angular  resolution  than  the  Discus  buoys,  but  it  sunk  prior  to  the 
ECONOMEX  deployments.  We  should  consider  here  the  possible  effect  of  using 
a  surface  buoy  moored  in  95  meters  of  water  (Discus  C)  on  the  predictions  of 
noise  at  the  450  meter  and  790  meter  OBS’s.  We  can  note  from  the  measured 
wave  spectra  that  there  is  little  energy  below  0.1  Hz  in  the  wave  spectrum.  If 
we  calculate  tanh(fcd)  for  the  worst  case  of  0.1  Hz  with  d  =  95  m  we  find  our 
dispersion  relationship  is  in  error  by  less  than  1%.  Thus  we  would  expect  no 
correction  need  be  made  to  our  wave  spectra  measured  at  Discus  C. 

Let  us  now  consider  the  angular  resolution  available  from  pitch  and  roll  mea¬ 
surements.  Longuet-Higgins  et  al.  (25]  were  the  first  to  investigate  the  angular 
response  of  these  buoys  as  follows.  If  the  wavelengths  of  the  surface  motion  are 
large  with  respect  to  the  buoy  diameter,  the  buoy  tends  to  have  the  same  motion 
and  orientation  as  the  surface.  Then  if  we  measure  the  vertical  displacement  and 
the  two  angles  of  pitch  and  roll,  we  will  have  three  time  series  which  represent 
the  vertical  displacement  f,  and  its  spatial  derivatives  d$/dxi  and  df/dxj.  Using 
the  notation  of  the  last  chapter,  we  may  represent  the  sea  surface  as  a  stochastic 
integral 

$•  =  9 1 J  exp(t(#c  •  x  —  at))dS.  (5.1) 

Then,  since  k  =  (/c  cos  7, /c  sin  7),  we  can  write  our  three  time  series  denoted  by 
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6,6,6  as 


6  =  9 if  exp(t'(#c  •  z  —  at))dS 
Js 

=  9f  /  *#c  cos  7  exp(t(#c  •  z  -  crt))dS 
Js 

£s  =  9?  [  t/c  sin 7  exp(t(#c  •  z  —  <rt))dS. 
^  Js 


(5.2) 


Next  we  cam  form  co-spectra  C{j{a)  and  quadrature  spectra  Qij {a )  from  the  tune 
series  6,  6 ,  and  we  find 


Cn{a)  =  Jq  F(o,'i)d') 1, 

C22(o)  =  /  k2  cos2  7-F(<x,7)  d'), 
CiS(a)  =  /c2sin27F(<r,Tf)d'Y, 
Cjs(<r)  =  /  /c2  cos  7  sin  7-F(o,  7)  c*7, 

rt* 

Q  «(<r)  =  /  /ccos7F(a,7)d7, 

jo 

f2r 

Qis{o)  =  Jr  #csin'YF(<r,7)d7, 


(5.3) 


where  F(cr,  7)  =  n(a)G(<r,  7)  is  the  frequency-directional  spectrum  of  the  surface 
elevation.  We  can  define  the  Fourier  coefficients  of  F(o, 7)  as 

1  r2*- 


«.(«,)  +  =  - /  V’F(a,i)  d7. 

7T  JO 


(5.4) 


where 

I  n=oo 

F{tr,  7)  =  -ao  +  51  (°»  cos(ri7)  +  &nsin(n7)].  (5.5) 

*  n=l 

We  can  see  that  the  right  hand  sides  of  Eq.  (5.3)  are  related  to  these  Fourier 
coefficients  as  follows: 

oq(o)  =  -Cn(o), 

JT 

ai{<r)  —  -^Qn{a)i  b2(o)  =  -^Q\i{cr), 

a2(a)  =  -^j~j\C22(o)  -  C'ss(a)],  b2(a)  =  -^j^C2z(o).  (5.6) 
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The  pitch  and  roll  buoy  then  gives  us  the  first  five  coefficients  in  the  Fourier  series 
describing  the  angular  spectrum  of  the  surface  elevation  at  each  of  the  frequencies 
for  which  we  find  the  co-spectra  and  quadrature  spectra  of  the  time  series.  We 
must  now  use  these  five  Fourier  coefficients  to  find  estimates  of  the  frequency- 
directional  spectra  ft  (a)  and  G(a,  7). 

Our  estimate  of  the  wave  power  spectrum 

ft(o)  =  f  F{p,7)  =  7roo(o)  =  Cn(a),  (5.7) 

Jo 

is  obvious  from  inspection  of  Eq.  (5.5),  as  well  as  our  definition  of  the  co-spectra. 

The  best  estimate  of  the  directional  spectrum  is  not  as  simple.  An  obvious 
choice  would  be  to  try  the  truncated  sum 

G(c r, 7)  =  H — —  (ai  cos  7  +  61  sin 7  +  02  cos 27  +  62  sin 27).  (5.8) 

27T  JTCIq 

This  sum  is  actually  a  convolution  of  the  true  directional  spectrum  with  a  weight¬ 
ing  function,  and  considerable  smoothing  results  in  the  estimated  spectrum.  Cal¬ 
culations  of  Ia33  made  using  this  type  of  directional  spectrum  estimate  from  simu¬ 
lated  directional  spectra  are  typically  in  error  by  a  factors  of  10s.  Other  weighted 
averages  of  the  first  five  Fourier  coefficients  can  be  made,  but  they  too  produce  a 
much  smoothed  estimate. 

Several  investigators  have  fit  empirical  curves  to  measured  directional  spectra. 
Longuet-Higgins  has  suggested  the  wave  directional  spectra  fit  the  form  [25] 

G(o,  7)  a  |  cos2'  (7/2)  |  (5.9) 

where  the  spreading  parameter  5  is  a  function  of  frequency  and  wind  speed.  Kib- 
blewhite  and  Wu  used  an  empirical  relationship  to  find  s  based  on  the  wind  speed 
and  frequency,  and  then  calculated  Ia33  analytically  [3].  One  could  also  match 
the  measured  first  five  Fourier  coefficients  to  the  first  five  Fourier  coefficients  of 
the  empirical  spectrum  cos2*  7/2  to  estimate  the  parameter  s ,  as  suggested  by 
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Longue t-Higgins  [25],  and  hence  calculate  Ias3-  However,  the  data  used  to  de¬ 
velop  the  empirical  formula  are  generally  taken  under  conditions  of  steady  wind 
speed  and  direction.  In  our  field  data,  the  wind  speed  and  direction  cam  vary 
significantly,  giving  rise  to  wave  fields  with  different  directionality,  and  in  gen¬ 
eral  a  broader  directional  spectrum  than  that  predicted  by  the  empirical  formula. 
Therefore,  we  would  expect  predictions  of  Ia33  calculated  from  spectra  derived 
from  empirical  formulas  to  be  lower  than  the  true  value  under  variable  meteoro¬ 
logical  conditions.  In  particular,  if  the  true  directional  spectrum  is  bimodal  (two 
peaks  corresponding  to  two  wave  fields  generated  by  winds  in  different  directions) 
there  can  be  significant  energy  in  opposing  wave  directions  G^G^  +  tt)  which  is 
not  predicted  by  the  cosine  power  curve.  Donelan  et  al.  have  suggested  a  better 
fit  to  the  data  is  found  in  a  sech2/3'y  distribution  (26]  with  /?  being  the  spreading 
parameter,  but  estimates  of  Ias3  based  on  estimates  of  (3  are  also  too  small. 

Another  approach  would  be  to  use  a  data  adaptive  spectral  estimation  tech¬ 
nique  such  as  the  maximum  likelihood  method  (MLM)  or  the  maximum  entropy 
method  (MEM).  The  method  used  here  to  estimate  the  directional  spectrum  given 
the  first  five  Fourier  coefficients  is  the  MEM.  This  method  produces  a  spectal  es¬ 
timate  which  retains  the  first  five  Fourier  coefficients  and  estimates  the  remaining 
coefficients  based  on  the  first  five.  We  will  follow  Lygre  and  Krogstad  [27]  in 
developing  an  algorithm  to  make  this  estimate. 

We  will  define  a  function  with  a  Fourier  series  (suppressing  the  dependence  on 
o)  on  the  interval  (— 7r,7r)  as 

D(l)  =  i  £  C\  CO  =  1,  c.„  =  <.  (5.10) 

n=-oo 

The  entropy  of  D  is  defined  by 

H(D)  =  1 log (D(7))^,  (5.11) 

and  it  has  been  shown  by  Burg  [28]  that  the  function  maximizing  H(D)  subject 
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to  the  constraint  that  the  coefficients  cn  equal  some  known  c*  for  k  <  N  is 


D(l)  = 


0\ 


2?r  |1  -  - <f>Ne~iN^D |2’ 

where  and  0*  are  obtained  from  the  Yule- Walker  equations 


(5.12) 


and 


1 

H 

l 

\? 

<i>  i 

Cl 

Cl 

•  • 

«  • 

•  • 

•  • 

c[ 

X 
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• 

cn- i  *  *  * 

Cl  1 

<t> N 

cn 

(5.13) 


o]  =  1  —  M - 


(5.14) 


In  our  case  we  have  N  =  2  with  Ci  =  (ai  +  t&i)/ao  and  c2  =  (a2  -f-t&2)/oo-  We  can 
now  solve  this  system  of  equations  to  find  <f>n  in  terms  of  cn: 


<f>  1  =  (ci  —  c2Cl)/(l  —  lcl|2)» 

<fa  =  «2  —  Cl^i,  (5.15) 


and  finally  we  can  substitute  these  into  Eq.  (5.12)  to  find  our  directional  spectral 
estimate 


D(  l)  = 


1  —  <f>xc\  —  <j>iC*2 
2jt|1  —  <f>ie~i'1  —  <f>2C2il\2' 


(5.16) 


We  now  have  an  estimate  of  the  angular  distribution  of  the  wave  energy  at  each 
frequency  which  reproduces  the  first  five  Fourier  coefficients  which  produced  it, 
and  uses  these  and  the  Yule- Walker  equations  to  extrapolate  the  remaining  co¬ 
efficients.  This  technique  has  been  shown  by  Lygre  and  Krogstad  [27]  to  give  a 
much  more  peaked  distribution  over  the  MLM  technique,  and  to  resolve  a  bimodal 
wavefield. 

Figure  3  shows  the  estimated  spectra  from  the  various  methods  for  a  simulated 
bimodal  spectral  input,  and  Fig.  4  shows  the  estimates  of  directional  spectra  for 
the  0.16  Hz  bin  from  Discus  E  on  January  27  at  1200.  This  datum  was  chosen  as 
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illistrative  because  it  occurs  after  a  shift  in  wind  direction,  and  one  could  expect 
the  true  spectrum  to  be  bimodal.  One  can  see  for  the  estimates  from  the  simulated 
data  that  the  MEM  estimate  provides  a  better  representation  of  the  structure  of 
the  actual  spectrum.  One  can  also  note  in  the  comparison  of  actual  data  a  bimodal 
structure  is  evident  in  the  MEM  estimate  only. 


5.2  Coupling  factors  and  bottom  gain 

While  the  coupling  factors  /  HfE^k  dk  can  be  evaluated  analytically  for  an  ocean 
of  infinite  radius,  [4]  they  predict  infinite  noise  when  the  horizontal  components 
are  considered.  We  can  solve  this  by  adding  a  relaxation  mechanism  and  thereby 
attenuation  in  the  water  column,  or  we  can  use  a  finite  effective  radius  with  some 
physical  basis  such  as  the  ocean  basin  radius,  storm  radius  (if  applicable),  or  such. 
The  coupling  factors  must  then  be  evaluated  either  approximately  analytically  or 
numerically.  The  attenuation  due  to  absorption  at  the  frequencies  of  interest  is 
believed  to  be  too  small  to  limit  the  noise  on  ocean  basin  scales,  and  so  we  will 
assume  an  effective  radius.  Calculations  show  that  at  distant  ranges  the  horizontal 
coupling  factors  vary  roughly  as  \/R,  and  thus  doubling  the  effective  radius  has 
only  a  minor  effect  on  the  overall  source  strength  Qo(w) 

Based  on  the  distance  to  shore  of  the  ECONOMEX  deployments  and  on  trans¬ 
mission  loss  studies  in  this  frequency  range,  an  effective  radius  of  100  kilometers 
was  chosen.  The  magnitude  of  the  coupling  factors  was  then  calculated  numeri¬ 
cally  for  receiver  depths  of  450  meters  and  2500  meters.  These  results  are  shown 
in  Fig.  5. 

The  final  step  in  making  our  predictions  is  the  calculation  of  the  bottom  gain, 
2£,  for  the  horizontal  and  vertical  components.  The  most  computationally  inten¬ 
sive  part  of  this  task  is  the  calculation  of  the  depth  dependent  Green’s  functions, 
g(k,  z,  z1),  for  which  we  use  the  fast  field  approach.  The  fast  field  approach  can  be 
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generalized  to  a  fully  elastic  elastic  media,  but  requires  that  the  wave  equation  be 
separable  in  depth,  that  is  to  say  range-independent.  We  will  use  the  fast  field  ap¬ 
proach,  that  realizing  errors  that  may  ensue  due  to  range  dependent  bathymetry 
and  media. 

The  tool  used  to  solve  the  full  wave  problem  in  this  work  is  the  SAFARI  set  of 
programs  (Seismo- Acoustic  Fast  field  Algorithm  for  Range  Independent  environ¬ 
ments)  developed  by  Schmidt  and  Jensen  [13,14].  The  basic  solution  technique  in 
SAFARI  is  to  represent  the  field  in  each  of  a  series  of  homogeneous  layers  by  the 
Hankel  transforms  of  the  unknown  potentials  satisfying  the  homogeneous  wave 
equation.  The  boundary  conditions  at  each  layer  interface  yield  a  set  of  local 
equations  involving  the  unknown  potentials  of  the  adjacent  layers.  These  local 
equations  are  collected  into  a  global  matrix  which  can  then  be  solved  to  yield  all 
the  unknown  potentials  simultaneously.  Solutions  are  determined  efficiently  by 
implementing  modern  numerical  techniques.  An  additional  advantage  relating  to 
this  work  is  the  inclusion  of  the  integration  of  Eq.  (4.4 1)  as  an  option  in  the  code. 

The  environmental  model  used  in  SAFARI  requires  a  number  of  parameters. 
The  user  must  specify  for  each  layer  the  compressional  wave  speed,  shear  wave 
speed,  compressional  and  shear  attenuations,  and  density.  The  uppermost  and 
lowermost  layers  are  taken  to  be  semi-infinite  half-spaces.  Due  to  the  long  wave¬ 
lengths  involved  at  the  very  low  frequencies  of  our  predictions,  the  environmental 
model  should  be  as  accurate  as  possible  fairly  deep  into  the  bottom.  It  is  fortu¬ 
nate  that  compressional  wave  speed  data  from  a  deep  seismic  reflection  study  of 
the  U.S.  mid-Atlantic  continental  margin  was  made  available  prior  to  publication. 
The  EDGE  seismic  experiment  [15]  involved  recording  seismic  profiles  to  16  sec¬ 
onds  off  the  Virginia  coast  in  the  same  region  as  the  ECONOMEX  and  SWADE 
experiments.  Using  the  compressional  wave  speed  data,  the  figures  and  equations 
of  Hamilton  [29]  and  discussions  with  other  investigators  working  in  the  area  of 
geoacoustic  modeling  [30,31],  geoacoustic  models  were  developed  for  the  bottom 
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layer 

depth  (m) 

Cp  (m/s) 

Ct  (m/s) 

1.  ( dB/X ) 

p  (s/cm3) 

vacuum 

* 

0 

0 

0 

0 

0 

fluid 

0 

1500 

0 

.015 

0 

1 

elastic 

2500 

1756 

200 

.5 

1 

1.8 

elastic 

3165 

2193 

540 

.3 

.5 

2.0 

elastic 

3678 

2495 

895 

.1 

.3 

2.2 

elastic 

5011 

2648 

1314 

.1 

.3 

2.4  | 

elastic 

5276 

3817 

1900 

.1 

.3 

2.7 

elastic 

7231 

4042 

2011 

.1 

- 

.3 

2.8 

elastic 

8896 

6380 

3180 

.05 

.1 

3.0 

Table  5.2:  Deep  site  environmental  model.  CPi ,  and  7Pi<  are  compressional  and 
shear  wave  speeds  and  attenuations,  respectively;  p  is  the  density. 


■US 

Eieai 

vacuum 

* 

0 

0 

0 

0 

0 

fluid 

0 

1500 

0 

.015 

0 

1 

elastic 

450 

1799 

200 

.5 

l 

1.8 

elastic 

970 

2174 

540 

.3 

.5 

2.0 

elastic 

1396 

2683 

895 

.1 

.3 

2.2 

elastic 

2441 

3441 

1720 

.1 

.3 

2.4 

elastic 

3332 

4398 

1739 

.1 

.3 

2.5 

elastic 

4938 

5872 

2933 

.1 

.3 

2.7 

elastic 

6582 

6159 

3059 

.05 

.1 

2.9 

elastic 

8227 

6380 

3170 

.05 

.1 

3.0 

Table  5.3:  Shallow  site  environmental  model.  CPi,  and  7Pi,  are  compressional  and 
shear  wave  speeds  and  attenuations,  respectively;  p  is  the  density. 


at  locations  corresponding  to  water  depths  of  450  meters  and  2500  meters.  The 
models  thus  developed  are  given  in  tables  5.2  and  5.3.  Depth  profiles  of  com¬ 
pressional  and  shear  wave  speeds  at  the  two  sites  are  given  in  Fig.  6.  SAFARI 
calculations  of  for  both  m  =  1  and  m  =  2  for  a  receiver  at  450  meters  and 
2500  meters  are  given  in  Fig.  7. 

We  can  obtain  a  qualitative  understanding  of  the  relative  importance  of  the 
different  propagation  mechanisms  by  looking  at  Fig.  8,  contour  plots  of  the  magni¬ 
tude  of  the  wavenumber  integrands  plotted  against  the  inverse  of  the  phase  speed, 
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or  “slowness”  k/oj,  and  frequency.  One  can  see  for  the  deep  site  the  normal  modes 
play  the  predominant  role  in  propagating  the  surface  noise  to  the  deep  receiver, 
as  only  for  the  lowest  frequencies  do  we  get  a  contribution  from  waves  with  phase 
speeds  less  than  1500  m/s.  In  the  shallow  case  there  is  a  significant  contribution 
at  all  the  frequencies  of  interest  from  the  lower  phase  speeds,  indicating  the  im¬ 
portance  of  interface  waves.  These  waves  are  propagating  horizontally  but  suffer 
an  exponential  decay  in  the  vertical,  and  thus  they  are  excited  only  in  the  shallow 
case. 

We  can  combine  the  effects  of  the  coupling  factors  and  the  bottom  gain,  along 
with  the  constants  in  Eq.  (4.44)  to  find  the  temporally  invariant  part  of  the 
solution,  and  we  can  write  our  prediction  as  the  one-sided  sound  pressure  spectrum 
level  (dB  re  1  nPa2/Hz) 

SL(f)  =  20  log  fl(//2)  +  lOlog /«„(/)  +  10  log  !?(/),  (5.17) 


where  f  is  the  acoustic  frequency  /  =  w/(2xj  and  B(f)  is  this  time  invariant  part: 


10log£(/)  =  197  +  10  log  / 


+ 


(z?(/)  /*■  % 

T’(f)  £' Bsn;kdk\\. 


/M*kdk 


(5.18) 


The  quantity  10  [og  B(f)  is  shown  in  Fig.  9  for  both  the  deep  and  shallow  sites. 
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Chapter  6 


Predictions  for  Receivers  at  450 
meters  and  2500  meters 

S.l  Variation  of  spectral  level  with  frequency 
(spectral  shape) 

Since  the  wave  height  power  spectrum  appears  in  Eq.  (4.19)  as  ft2(/ /2),  we  might 
expect  the  shape  of  the  acoustic  spectrum  to  be  related  to  the  shape  of  the  wave 
height  spectrum  at  double  the  wave  frequency,  but  we  will  see  this  shape  is  mod¬ 
ified  by  a  number  of  factors.  The  wave  height  spectrum  generally  shows  a  very 
steep  rise  to  a  spectral  peak  followed  by  a  somewhat  gentler  (a  a-4)  slope  at 
frequencies  above  the  peak  [26].  This  peak  frequency  is  generally  characterized  as 
being  inversely  proportional  to  wind  speed,  and  thus  is  usually  lower  in  frequency 
at  higher  wave  heights.  The  peak  tends  to  be  quite  narrow,  normally  occurring 
at  frequencies  of  about  0.1  to  0.2  Hz,  although  in  a  newly  developing  wave  field 
it  can  be  higher.  This  tends  to  give  rise  to  an  acoustic  spectrum  with  a  peak  in 
the  0.2  to  0.4  Hz  range. 

This  overall  shape  will  of  course  be  modified  by  the  effects  of  the  directional 
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spectrum,  the  coupling  factors,  and  the  bottom  gain  as  functions  of  frequency. 
Effects  of  the  coupling  factors  are  easily  seen  in  Fig.  5  and  tend  to  emphasize  the 
higher  frequencies,  which  will  tend  to  mitigate  the  slope  of  the  peak  in  the  acoustic 
spectrum  at  frequencies  above  the  peak.  The  spectral  slope  of  the  coupling  factors 
is  oc  oj2.  The  bottom  contribution  to  the  spectral  shape  is  more  complicated,  as 
seen  in  Fig.  7.  In  the  shallow  case,  a  peak  in  the  bottom  gain  appears  at  about 
0.18  Hz,  which  sharply  drops  to  a  low  at  about  0.28  Hz,  and  follows  with  a  rise 
above  0.28  Hz  with  a  slope  of  about  u2.  This  will  tend  to  flatten  the  acoustic 
spectrum  if  the  peak  in  the  wave  spectrum  occurs  above  0.1  Hz  (as  it  almost 
always  does).  The  peak  in  the  bottom  gain  at  0.18  Hz  is  normally  well  overcome 
by  the  sharp  drop  in  the  wave  height  spectrum  at  wave  frequencies  below  0.1  Hz. 
In  the  deep  case,  the  peak  in  the  bottom  gain  occurs  at  about  0.24  Hz,  and  is  more 
likely  to  have  an  effect  on  the  shape  of  the  acoustic  spectrum.  In  most  cases,  the 
frequency  of  this  peak  in  the  bottom  gain  falls  slightly  below  twice  the  frequency 
in  the  wave  height  spectrum,  thus  serving  to  broaden  the  peak  in  the  acoustic 
spectrum.  Above  0.3  Hz  the  bottom  gain  is  relatively  fiat  in  the  deep  case.  Thus 
in  the  deep  location  the  bottom  will  tend  to  enhance  the  peak  in  the  acoustic 
spectrum,  and  at  frequencies  above  the  peak  the  down  slope  should  be  greater  at 
the  deep  location  than  in  the  shallower  case.  When  we  look  at  the  overall  transfer 
function  B(f)  shown  in  Fig.  9,  we  see  for  the  shallow  site  a  post  peak  slope  of 
oc  w3  and  for  the  shallow  case  a  slope  of  oc  w5. 

By  far  the  greatest  effect  on  the  shape  of  the  noise  spectrum  next  to  the  wave 
height  power  spectrum  is  the  behavior  of  /*„.  If  one  examines  the  behavior  of 
the  empirical  models,  it  is  clear  that  the  value  of  the  spreading  parameters  s  or 
P  vary  such  as  to  reduce  the  value  of  Jass  at  the  frequency  of  the  spectral  peak. 
Mitsuyasu  tt  al.  [32]  have  proposed  a  model  for  frequencies  above  the  spectral 
peak  in  which  s  varies  as  s  a  {U/ga)'26,  where  U  is  the  wind  speed.  In  Donelan’s 
model,  P  is  dependent  on  o/ap  only,  with  ap  being  the  frequency  at  the  peak  [26]. 
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Variation  of  Ia33  versus  frequency  for  this  model  is  shown  in  Fig.  10.  One  can  see 
that  for  either  of  these  empirical  models  the  directional  spectrum  is  narrower  near 
the  peak  frequencies,  and  broader  away  from  the  peaks,  and  thus  the  spreading 
integral  Ia33  will  be  a  minimum  near  the  peak.  The  effects  of  the  wave  height  power 
spectrum  and  the  spreading  integral  will  oppose  one  another,  and  again  tend  to 
flatten  the  spectrum.  We  can  see  in  Figs.  11-  26  the  estimates  of  10  logP(/), 
fl(a/2x),  20log(n(//2)),  and  10 log  /<*«(/)  for  several  illustrative  examples  for  a 
2500  m  receiver,  and  in  Figs.  27  -  30  the  same  quantities  for  a  450  m  receiver.  One 
can  see  from  Figs.  11-  30  that  the  wave  height  power  input  can  vary  by  as  much  as 
60  dB  in  a  given  spectrum  over  the  frequency  range,  whereas  the  spreading  input 
varies  over  a  much  narrower  range  of  up  to  20  dB  in  a  given  spectrum.  This  will 
cause  the  wave  height  power  spectrum  to  dominate,  giving  rise  to  a  spectral  peak, 
although  one  much  diminished  from  the  peak  in  the  generating  wave  spectrum. 

In  summary  we  find  the  shape  of  the  acoustic  spectrum  will  in  general  resemble 
that  of  the  generating  wave  spectrum  in  that  their  will  be  a  sharp  rise  to  a  spectral 
peak,  but  in  the  acoustic  spectrum  the  peak  will  be  of  a  lower  magnitude  and 
broader  than  that  present  in  the  wave  spectrum. 


6.2  Variation  of  spectral  level  with  wave  height 

The  variation  in  the  predicted  noise  with  wave  height  will  be  the  result  of  the 
combined  effects  of  the  wave  height  spectrum  squared  and  the  spreading  integral. 
We  saw  in  the  last  section  that  the  two  effects  opposed  one  another,  with  the 
variation  in  the  wave  height  spectrum  dominating,  giving  rise  to  a  spectral  peak. 
When  we  look  at  predicted  spectra  from  different  times  corresponding  to  different 
meteorological  conditions,  we  find  the  variation  in  the  wave  height  spectra  is  again 
greater  than  that  of  the  spreading  integral.  We  would  thus  predict  an  increase  in 
noise  level  with  an  increase  in  overall  wave  height  or  sea  state,  with  the  increase 
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in  wave  height  again  partially  offset  by  the  decrease  in  the  spreading. 

To  allow  us  to  see  the  effects  of  increasing  wave  height  on  our  predictions, 
it  is  best  to  study  a  period  of  relatively  constant  wind  direction  with  increasing 
wind  speed  and  corresponding  wave  height.  This  will  minimize  the  effects  of 
other  conditions  which  could  effect  the  wave  directional  spectrum,  such  as  wind 
direction  versus  fetch  direction,  sharp  changes  in  wind  direction,  etc.  We  can  see 
from  Fig.  31  that  we  have  such  a  time  available  to  us  at  Discus  E  from  0300  to 
1200  on  January  21,  1991.  During  that  time,  wind  speed  rises  from  about  5  m/s 
to  about  12  m/s  with  a  rise  in  significant  wave  height  of  0.25  m/s  to  1.9  m/s. 
Some  predicted  noise  spectra  for  this  period  are  included  as  Figs.  11-14.  One  can 
see  in  the  noise  prediction  at  0300  a  peak  forming  at  0.6  Hz  with  a  level  of  121  dB. 
In  subsequent  predictions,  this  peak  moves  lower  in  frequency,  to  0.45  Hz  at  0600, 
0.38  Hz  at  0900  and  finally  0.30  Hz  at  1200.  The  peak  level  also  increases  from 
121  dB  to  125  dB,  130  dB,  and  finally  133  dB  at  1200.  It  is  interesting  to  note 
the  wave  height  power  at  the  peak  for  this  period  varies  from  -13  dB  to  12  dB,  a 
range  of  25  dB,  while  the  value  of  the  spreading  integral  at  the  peak  varies  from 
-18  dB  to  -25  dB,  for  a  range  of  -7  dB.  This  clearly  shows  the  mitigating  effect 
of  the  spreading  integral  in  limiting  the  noise  power  at  the  peak  under  moderate 
conditions. 

It  is  also  interesting  to  consider  the  predictions  in  the  case  of  severe  weather 
conditions.  The  highest  significant  wave  height  recorded  for  which  ECONOMEX 
data  is  available  was  5.9  meters  at  Discus  E  in  the  late  afternoon  of  March  4. 
Again  the  wind  direction  was  relatively  steady.  The  meteorological  data  for  this 
period  is  shown  in  Fig.  32,  and  the  noise  predictions  for  0400  to  1600  on  that 
day  are  seen  as  Figs.  15-21.  The  peak  noise  level  prediction  for  this  time  is  155 
dB,  with  a  wave  height  input  of  35  dB  and  a  spreading  integral  value  of  -21  dB. 
Investigation  of  our  MEM  estimates  of  shows  that  it  reaches  a  minimum  of 
about  -25  dB  at  the  peak  frequency  under  conditions  of  moderate  wave  height  (~  3 
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meters),  with  no  further  decrease  with  increasing  wave  height  above  that  level. 
This  means  our  estimated  directional  distribution  is  not  becoming  increasingly 
narrow  at  wave  heights  above  ~  3m.  This  is  in  disagreement  with  the  observed 
behavior  of  the  directional  spectra  of  wave  height  at  high  sea  state  (26).  We 
might,  therefore,  expect  our  predictions  to  overestimate  the  noise  at  high  sea 
states.  The  increase  in  predicted  noise  level  at  wave  heights  above  this  threshold 
will  be  directly  proportional  to  the  increase  in  the  wave  height  power  squared. 

At  frequencies  above  the  peak,  the  variation  in  wave  height  power  with  wave 
height  is  much  less  than  that  at  the  peak.  One  can  see  that  the  variation  of  the 
average  wave  height  power  on  January  27  from  0300  to  1200  in  the  0.5  Hz  to 
0.68  Hz  was  only  from  ~  -20  dB  to  ~  -15  dB.  This  is  typical  of  most  of  the  data 
analyzed  in  this  frequency  range.  Even  under  the  extreme  conditions  of  March  4 
the  wave  height  power  in  this  frequency  range  does  not  increase  much  above  this 
level.  Similarly  the  spreading  input  varies  little  in  this  range  of  frequencies  for  a 
steady  wind  direction,  varying  from  — 5  dB  to  ~  0  dB.  We  would  therefore  predict 
very  little  variation  in  the  noise  level  in  this  frequency  range  for  an  increasing  sea 
state. 


6.3  Variation  of  spectral  level  with  changing  wind 
direction 

One  of  the  more  interesting  studies  we  can  make  involves  examining  the  varia¬ 
tion  of  noise  level  under  conditions  of  constant  wind  speed  but  changing  wind 
direction.  A  theory  which  proposes  that  wind  turbulence  is  the  direct  cause  of 
the  acoustic  noise  would  predict  very  little  variation  in  noise  levels  under  these 
conditions.  Thus  the  predictions  we  make  here  may  be  useful  in  judging  the  ac¬ 
tual  contribution  from  the  different  mechanisms.  Under  our  theory,  we  intuitively 
expect  the  noise  to  increase  from  the  increased  spread  in  the  directional  spectrum 
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as  the  new  wave  field  is  developed,  but  we  also  realize  the  original  wave  field  is 
diminishing,  due  to  the  loss  of  wind  forcing  and  the  combined  effects  of  dissipation 
and  non-linear  interaction^  with  «,he  newly  developing  wave  field. 

We  can  see  these  combined  effects  in  the  data  from  Discus  E  in  the  early  hours 
of  February  23,  1991.  As  we  can  see  in  Fig.  33  the  wind  direction  veers  sharply 
by  about  60  degrees  at  0200  while  the  wind  speed  stays  relatively  constant  at  8 
to  11  m/s.  The  response  in  the  predicted  noise  spectra  is  shown  in  Figs.  22-26. 
One  can  see  the  noise  prediction  at  the  peak  is  fairly  constant  at  about  135  dB 
from  0200  through  0400,  then  begins  to  rise  until  it  reaches  140  dB  at  0600.  The 
wave  height  power  input  during  this  period  actually  drops  from  15  dB  to  8  dB 
at  0400,  and  it  rises  back  to  about  15  dB  at  0600.  The  spreading  input  at  the 
peak  increases  from  about  -23  dB  at  0200  to  -17  dB  at  0400,  and  then  remains 
at  about  this  value  through  0600.  Thus  the  increase  in  predicted  noise  is  due  to 
the  increased  spread  in  the  directional  spectrum,  not  an  increase  in  wave  height. 
Another  effect  of  the  wind  shift  on  the  predicted  acoustic  spectrum  is  that  the 
peak  tends  to  broaden.  This  is  seen  most  clearly  at  0400  in  Fig.  24,  when  the 
peak  in  the  developing  wave  field  and  the  peak  in  the  pre-existing  wave  field  are 
both  of  similar  magnitudes  but  different  frequencies.  This  behavior  is  typical  of 
the  wind  shifts  analyzed. 


6.4  Variation  of  spectral  level  with  receiver  depth 

The  variation  in  predicted  spectral  levels  with  receiver  depth  can  be  attributed 
to  three  factors.  First,  for  a  receiver  in  very  shallow  water,  the  gravity  wave 
dispersion  relationship  will  begin  to  depart  from  Eq.  (4.4)  and  the  bottom  will 
begin  to  have  an  effect  on  the  wave  height  spectrum.  For  our  receivers  at  2500  m 
and  450  m,  this  is  not  the  case.  The  two  remaining  factors,  the  coupling  factor 
differences  and  the  bottom  gain  differences,  will  cause  variation  in  our  expected 
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noise  levels  at  these  two  depths. 

Referring  to  Fig.  5,  we  can  see  that  the  difference  in  the  coupling  factor  gain 
between  the  two  locations  varies  from  about  15  dB  at  0.06  h z  to  about  7  dB  at  0.68 
Hz,  with  the  shallow  site  having  the  higher  value.  Tne  bottom  gain  is  again  more 
complicated,  with  the  shallow  bottom  gain  being  higher  at  frequencies  less  than 
0.2  Hz  and  greater  than  0.5  Hz,  and  the  deep  bottom  gain  being  higher  between 
about  0.2  Hz  and  0.4  Hz.  The  combined  effects  are  visible  in  Fig.  9,  where  we  can 
see  the  differences  in  the  energetic  part  of  the  spectrum  are  quite  minimal.  Given 
the  same  wave  height  directional  spectrum  input,  we  expect  the  deep  case  to  yield 
higher  noise  levels  of  about  6  dB  at  0.25  Hz,  and  we  expect  the  shallow  case  to 
yield  higher  noise  levels  of  about  5  dB  at  0.68  Hz.  Due  to  the  limited  amount  of 
data  analyzed  and  the  variance  of  the  estimated  spectra,  it  is  difficult  to  see  this 
variation  in  the  predicted  spectra  between  the  two  depths.  We  will  therefore  not 
try  to  judge  the  success  of  the  predictions  in  this  area. 
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Chapter  7 

Experimental  Noise 
Measurements 


7.1  Experimental  description 

ECONOMEX  (Environmentally  Controlled  Oceanfloor  Noise  Monitoring  Experi¬ 
ment)  [10]  was  designed  to  provide  a  long-term,  high  quality  seismo-acoustic  noise 
data  set  which  could  be  coupled  to  the  surface  wave  and  meteorological  data  of 
the  SWADE  experiment.  The  instrumentation  consisted  of  six  Office  of  Naval 
Research  (ONR)  ocean  bottom  seismometers  (OBS’s)  and  two,  one  vertical  and 
one  horizontal,  75  meter  six  element  hydrophone  arrays.  The  instruments  were  off 
the  Virginia  coast  in  January  1991,  recovered  in  February  1991  for  maintenance 
,  and  redeployed  from  February  through  early  April  1991.  Precise  instrument 
locations  and  deployment  dates  are  listed  in  table  7.1. 

The  ONR  OBS  instruments  deployed  consisted  of  a  three-component  geophone 
for  measuring  ground  motion  in  the  0.07  to  80  Hz  range,  a  Cox- Webb  differential 
pressure  gauge  (DPG)  for  measuring  long  period  pressure  signals  in  the  water 
column,  and  in  the  original  deployment,  an  OAS  hydrophone  [34].  In  the  later 
deployment,  these  hydrophones  were  removed  to  improve  instrument  reliability. 
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Instrument 

Lat.  (N) 

Lon.  (W) 

Depth 

Deployment  dates 

Vert,  array 

2573m 

Jan  25-Feb  6;  Feb  22- Apr  5 

Hor.  array 

37°  24.7' 

Kmsmm 

2573m 

Jan  25-Feb  6;  Feb  22-Apr  5 

OBS  56 

37°  24.7' 

mkmxm 

2548m 

Jan  25-Feb  6;  Feb  22-Apr  5 

OBS  58 

37°  26.4' 

73°  31.4' 

2417m 

Jan  24-Feb  7;  Feb  22-Apr  5 

OBS  61 

*W**I31 

73°  24.4' 

2600m 

Jan  25-Feb  6;  Feb  21-Apr  4 

lITajzKFsy^ 

37°  33.2' 

mgamsm 

769m 

Jan  11-Feb  6;  Feb  22-Apr  5 

OBS  63 

37°  54.1' 

7 4°  16.5' 

443m 

Jan  10-Feb  6;  Feb  22-Apr  5 

OBS  51 

37°  35.$ 

74°  213' 

95m 

Jan  10-Feb  6;  Feb  22-Apr  5 

Table  7.1:  ECONOMEX  instrument  summary.  Positions  and  depths  listed  are  for 
the  second  leg  of  the  experiment;  those  for  the  first  leg  differ  only  slightly.  Also 
note  OBS  frames  61  and  62  exchanged  positions  between  the  first  and  second  legs, 
although  the  instruments  on  them  were  exchanged  also  such  that  the  instruments 
remained  deployed  in  the  same  locations. 


The  sensors  were  connected  via  pi  .  amplifiers  to  an  acquisition  package  consisting 
of  a  pre-whitening  and  anti-aliasing  filter,  a  gain-ranging  amplifier  to  improve 
dynamic  range,  and  an  analog  to  digital  converter.  The  combined  filter  response 
is  shown  in  Fig.  34.  The  acquisition  package  fed  a  recording  package  consisting  of  a 
.  RAM  buffer  and  an  optical  disc  recording  system  capable  of  storing  400  megabytes 
of  data.  The  typical  OBS’s  were  programmed  for  continuous  8  Hz  recording,  with 
the  anti-aliasing  filter  set  to  2  Hz;  however  OBS  56  was  set  to  record  at  128  Hz 
with  its  anti-aliasing  filter  set  to  40  Hz. 

The  75  meter  horizontal  and  vertical  arrays  each  consisted  of  six  OAS  hy¬ 
drophones  at  15  meter  separation.  The  hydrophone  signals  were  preamplified  by 
a  low-noise,  wide-range  preamplifier  and  sent  to  acquisition  and  recording  pack¬ 
ages  identical  to  those  of  the  OBS’s,  with  128  Hz  sampling  and  the  anti-aliasing 
filter  set  to  40  Hz.  In  the  second  deployment  the  bottom  three  hydrophones  of 
the  vertical  array  were  not  included  due  to  a  cable  malfunction.  The  array  cable 
jacket  included  loose  ended  fiber  strands  to  reduce  strumming  noise. 
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7.2  Data  selection  and  processing 


Following  instrument  recovery,  the  ECONOMEX  data  were  transcribed  from  the 
optical  disc  to  magnetic  tape.  While  the  data  has  not  yet  been  transcribed  from 
the  binary  machine  format  to  a  standard  format  for  further  dissemination,  it 
is  possible  to  read  the  binary  format  and  produce  ASCII  files  for  limited  time 
periods.  The  present  work  has  concentrated  on  using  the  differential  pressure 
gauge  (DPG)  data  at  the  8  Hz  sampling  rate  to  minimize  the  data  processing 
involved  while  still  adequately  sampling  the  frequency  band  of  interest.  The  DPG 
data  were  used  instead  of  the  geophone  data  because  of  current  uncertainties  in 
the  geophone  response.  Since  the  model  developed  is  only  valid  for  receivers  at 
depths  greater  than  400  meters,  it  was  decided  to  analyze  selected  data  from  one 
deep  DPG  (OBS  58  at  2417  m)  and  the  450  m  DPG  (OBS  63).  This  allowed  the 
maximum  depth  variation  comparison  given  the  instrument  deployment  depths. 

Estimated  spectra  were  generated  by  removing  the  mean  and  any  linear  trend 
from  64  second  segments  of  the  time  series,  and  then  averaging  512  poim  fast 
Fourier  transforms  of  34  minute  sections  of  data  using  a  Hanning  window.  The 
response  of  the  pre-whitening  and  anti-aliasing  filter  was  then  removed.  The 
bandwidth  of  64  seconds  results  in  a  frequency  resolution  of  0.0156  Hz,  which 
is  consistent  with  the  spectral  resolution  of  our  predictions,  which  is  0.02  Hz. 
An  example  of  a  full  spectral  range  observation  in  shown  in  Fig.  35.  Since  the 
predictions  are  limited  to  a  0.06  to  0.68  Hz  band,  further  observed  spectra  shown 
in  this  work  are  limited  to  the  same  band  for  clarity.  For  the  most  part,  data  were 
analyzed  which  corresponded  to  the  same  time  periods  for  which  predictions  were 
made  in  Chap.  6. 
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7.3  Observed  results  and  comparison  with  pre¬ 


dictions 

7.3.1  Spectral  shape  and  overall  spectral  noise  levels 

In  general,  there  is  good  overall  agreement  between  the  predicted  spectral  levels 
at  most  frequencies  in  the  band  of  interest.  Fig.  36  shows  an  example  of  one  such 
case  at  the  2500  meter  site,  and  Fig.  37  shows  similar  results  at  the  450  meter  site. 
While  the  levels  at  a  given  frequency  may  differ  between  predicted  and  observed 
by  up  to  5  dB,  there  is  an  overall  correspondence  between  the  two.  It  is  interesting 
to  note  in  many  of  the  comparisons  between  observed  and  predicted  spectra  that 
small  peaks  exist  in  the  observed  levels  which  are  present  in  the  predictions  at 
the  same  frequenies  but  with  diffemt  magnitudes.  These  two  figures  represent 
examples  of  the  best  agreement  between  predicted  and  observed  spectra. 

More  typical  of  the  level  of  agreement  are  Figs.  38-43,  where  we  can  see  close 
agreement  at  frequencies  around  the  spectral  peak,  but  differences  away  from  the 
peak  of  up  to  7  dB  in  the  deep  case,  and  up  to  17  dB  in  the  shallow  case.  Here 
the  overall  shape  is  correctly  predicted,  but  the  peak  is  broader  or  narrower  in 
the  observed  spectra,  giving  rise  to  large  differences  in  the  high  slope  region.  In 
the  predictions  from  the  shallow  site,  there  is  a  tendency  to  predict  levels  that  are 
too  high  in  the  band  0.1  Hz  to  0.2  Hz.  Since  this  band  corresponds  to  the  peak  in 
the  bottom  response  in  the  shallow  case,  errors  in  the  bottom  model  may  account 
for  this  difference. 

7.3.2  Variation  with  wave  height 

The  observations  show  an  increase  in  noise  level  with  wave  height,  with  good 
agreement  between  observations  and  predictions  at  moderate  (~l-3  meters)  wave 
heights,  as  seen  in  Figs.  39-43.  At  very  low  predicted  noise  levels,  which  corre- 
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spond  to  times  of  low  wave  height,  there  is  a  disparity  between  predictions  and 
observations.  An  example  of  this  is  the  developing  wave  field  of  the  early  hours 
of  January  27  at  Discus  E.  As  one  can  see  in  Figs.  44-  46,  the  agreement  is  good 
at  frequencies  corresponding  to  the  peak  in  the  predictions  as  the  wave  field  de¬ 
velops,  but  the  low  levels  predicted  away  from  the  peak  are  not  confirmed  by 
observation.  By  1200  the  wave  field  is  developed,  and  Fig.  39  shows  the  good 
agreement  between  the  prediction  and  observation  at  this  time.  A  possible  cause 
of  the  error  in  the  predictions  at  low  wave  heights  is  the  assumption  of  spatial  ho¬ 
mogeneity  of  the  source  wave  field.  Under  very  low  local  wave  height  conditions, 
it  is  possible  for  a  much  stronger  wave  field  at  some  distance  to  dominate  the  noise 
field,  thus  making  the  predictions  made  from  the  local  wave  field  very  much  in 
error.  Another  possible  cause  of  the  differences  under  low  wave  height  conditions 
is  the  existence  of  another  source  mechanism  generating  acoustic  energy,  whose 
noise  is  normally  dominated  by  that  caused  by  orbital  motion.  Under  low  source 
strength  conditions  for  the  orbital  motion  noise,  this  assumed  source  may  now 
dominate,  giving  rise  to  the  errors  noted  above.  The  predictions  in  general  are  in 
reasonably  close  agreement  when  measured  significant  wave  height  is  above  ~  1 
meter.  For  the  period  of  the  ECONOMEX  data,  roughly  77  percent  of  the  wave 
height  measurements  are  above  this  threshold. 

We  can  see  from  the  data  of  March  4,  as  seen  in  Fig.  47-50,  the  predictions  again 
begin  to  deviate  at  very  high  wave  height  conditions  which  correspond  to  strong 
winds  at  relatively  constant  direction,  with  the  predictions  being  higher  than  the 
observations.  The  disagreement  at  frequencies  corresponding  to  the  peak  is  up  to 
10  dB  under  these  conditions.  There  is,  however,  the  same  general  trend  in  the 
observed  noise  data  as  exists  in  the  predictions,  that  of  higher  levels  at  higher 
wave  heights  and  wind  speeds.  This  would  tend  to  indicate  the  wave  directional 
spectral  estimate  provided  by  the  MEM  technique  is  overestimating  the  spread  at 
these  high  wave  heights.  The  empirical  models,  which  are  based  on  strong,  steady 
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winds,  predict  less  spreading  under  these  conditions  than  do  the  MEM  estimates. 
Thus,  under  these  conditions,  the  empirical  model  estimate  for  the  spreading  may 
be  more  accurate.  The  observations  match  the  predictions  fairly  well  in  the  region 
of  frequencies  above  the  peak  under  high  wave  height  conditions. 

7.3.3  Variation  with  changes  in  wind  direction 

The  predictions  agree  quite  well  with  the  observed  data  under  conditions  of  chang¬ 
ing  wind  direction.  We  can  again  consider  the  wind  shift  of  February  23  at  0200, 
depicted  in  Fig.  33,  as  typical.  The  observed  and  predicted  noise  levels  for  sub¬ 
sequent  times  are  shown  in  Figs.  51-  55.  The  observed  noise  spectrum  changes 
with  time  roughly  as  the  predicted  spectrum  does.  Of  particular  note  is  the  broad 
peak  in  the  observed  spectrum  at  0400,  and  the  higher  levels  seen  at  the  peak 
from  0500-0700. 

It  is  also  interesting  to  note  the  predictions  and  observations  during  a  time 
of  highly  variable  conditions  such  as  those  of  the  afternoon  of  January  27,  again 
depicted  in  Fig.  31.  The  observations  and  predictions  are  seen  in  Figs.  36,56,57. 
Again  we  see  relatively  close  agreement  between  predictions  and  observations, 
with  the  increased  spreading  adding  to  the  noise  as  the  wave  height  diminishes  to 
keep  the  overall  noise  level  fairly  constant. 
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Chapter  8 


Conclusions  and 
Recommendations  for  Future 
Research 


8.1  Conclusions 

The  most  important  conclusion  one  can  draw  from  this  work  is  the  apparent  im¬ 
portance  of  the  orbital  motion  contribution  to  the  total  noise  level  in  the  frequency 
band  0.1  to  0.7  Hz.  The  close  correspondence  between  the  predicted  and  observed 
spectra  under  a  wide  range  of  conditions  is  strong  evidence  that  the  true  noise 
generation  mechanisms  in  this  band  sire  dominated  by  the  interaction  of  opposing 
surface  gravity  waves.  The  fact  that  the  predictions  hold  under  the  conditions  of 
changing  wind  direction  but  constant  wind  speed  would  tend  to  negate  the  impor¬ 
tance  of  the  direct  input  of  the  wind  turbulence  on  the  sea  surface  as  an  important 
sound  generation  mechanism  under  most  conditions.  There  is  a  possibility  that 
some  other  source  mechanism  or  mechanisms  contribute  significantly  to  the  noise 
spectrum  in  this  band  under  conditions  of  low  wave  height. 

The  second  conclusion  one  may  draw  from  this  work  is  the  importance  of  the 
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measured  directional  spectrum  in  determining  the  overall  noise  level.  Any  model 
based  on  empirical  relationships  using  wind  speed  as  their  input  would  have  to  also 
include  the  growth  and  decay  of  the  wave  field  based  on  changing  wind  direction. 
The  models  which  predict  directional  wave  spectra  from  a  knowledge  of  wind  speed 
alone  would  be  doomed  to  failure  in  the  general  case  due  to  their  assumption  of 
a  fully  developed  wave  field  and  the  changing  nature  of  the  true  wave  field. 

Another  important  point  brought  out  by  this  work  is  the  importance  of  the 
propagation  mechanisms  in  determining  the  overall  levels  and  shape  of  the  acous¬ 
tic  noise  spectrum  in  this  frequency  band,  which  is  consistant  with  the  work 
of  Schmidt  and  Kuperman  [23].  It  is  imperative  that  one  take  into  account  the 
propagation  if  one  wishes  to  compare  source  levels  between  two  different  locations. 
Otherwise  differences  in  the  bottom  contribution  could  cloud  important  correla¬ 
tions  or  lead  to  incorrect  conclusions.  A  key  element  in  research  of  this  type  is 
the  geoacoustic  bottom  model,  and  the  availability  of  measured  wave  speeds  deep 
into  the  bottom  in  the  present  work  was  extremely  fortunate. 

8.2  Recommendations  for  future  research 

As  the  predictions  of  acoustic  noise  depend  so  critically  on  the  angular  spread  of 
the  wave  energy,  a  fruitful  line  of  investigation  would  be  to  obtain  directional  wave 
height  spectra  of  greater  angular  resolution  in  future  experiments.  It  is  unfortu¬ 
nate  that  the  Spar  buoy,  the  sensor  with  the  greatest  directional  capability  in  the 
SWADE  project,  was  lost  prior  to  the  ECONOMEX  deployments.  The  increased 
directional  resolution  could  have  been  used  not  only  in  estimating  spectra  at  the 
location  of  the  buoy,  but  also  could  have  helped  in  determining  the  error  in  the 
estimate  of  the  directional  spectra  made  when  we  use  the  MEM  technique  at  other 
locations.  In  lieu  of  higher  resolution  directional  wave  spectra,  an  attempt  could 
be  made  to  use  some  of  the  other  SWADE  data,  such  as  directional  spectral  es- 
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timates  made  by  shipborne  arrays  or  radar  backscatter,  to  improve  the  estimates 
from  the  pitch  and  roll  buoys. 

Another  investigation  which  should  provide  interesting  results  would  be  the 
correlation  of  large  amounts  of  the  ECONOMEX  data  with  the  various  SWADE 
parameters  such  as  wind  speed  and  wave  height,  and  perhaps  to  automate  the  gen¬ 
eration  of  acoustic  predictions  from  the  SWADE  directional  spectra.  The  present 
research  only  scratched  the  surface  of  the  data  available  from  this  long  term  ex¬ 
periment.  It  would  be  useful  to  find  the  amount  of  long  term  agreement,  and  even 
more  interesting  to  find  other  periods  of  disagreement  between  the  predictions 
made  with  the  model  developed  herein  and  ECONOMEX  observations. 

The  final  recommendation  for  research  in  this  area  would  be  the  investigation 
of  the  range  dependent  aspects  of  the  problem.  There  is  range  dependence  in 
the  both  the  source  mechanism  and  the  propagation  mechanisms,  which  perhaps 
gives  rise  to  some  of  the  disagreement  between  our  predictions  and  observations. 
The  shallow  site,  in  particular,  would  benefit  from  a  consideration  of  its  range 
dependent  bathymetry.  Buckingham  [33]  has  found  the  noise  in  a  wedged-shaped 
ocean  with  pressure-release  boundaries  to  closely  approximate  the  noise  field  in 
the  range  independent  case,  but  further  work  is  needed  to  extend  his  work  to  a 
more  general  bottom.  As  well,  the  investigation  of  the  spatial  variation  of  the 
wave  height  spectrum  should  prove  feasible  once  the  entire  SWADE  project  data 
is  collected  and  correlated.  This  could  then  be  used  to  more  accurately  predict 
the  noise  at  the  ECONOMEX  sites. 
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ECONOMEX  and  SWADE  instrument  locations 


Longitude 


Figure  1:  Approximate  locations  of  the  ECONOMEX  and  relevant  SWADE  instru¬ 
ments.  denotes  pitch  and  roll  buoy  location,  *x”  denotes  location  of 
one  ONR  OBS  and  two  hydrophone  arrays,  and  *o*  denotes  the  location 
of  one  ONR  OBS. 
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Figure  2:  Schematic  arrangement  of  ECONOMEX  and  SWADE  deployments. 
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Magnitude  of  G  Magnitude  of  G 


Figure  3:  Estimates  of  the  directional  spectrum  of  G( 7)  =  AT(cos8(Tf/2)+cos10('y/2+ 
ir/4)),  where  7  is  the  asimuthal  angle,  and  N  is  chosen  to  normalise  the 
spectrum.  The  MEM  estimate  was  made  using  the  Lygre-Krogstad  al¬ 
gorithm  mentioned  in  the  text.  The  MLM  estimate  was  made  using  an 
algorithm  by  Lacoss  [35 j.  The  LCS  estimate  refers  to  the  weighted  aver¬ 
age  of  the  first  five  Fourier  coefficients  suggested  by  Longuet-Higgins  ef 
al.  [25]  The  empirical  estimate  was  made  by  estimating  the  parameter 
a  from  the  first  five  Fourier  coefficients  and  using  the  empirical  formula 
G(  7)  =  iV  cos3*  (7/2). 
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Directional  spectral  estimates  for  Discus  E  at  0.16  Hz,  1/27/91  1200 


Coupling  factor  integrals;  2500  m  and  450  m  receivers 


Figure  5:  Plots  of  10 log J*  fall;  +  H2H^)kdk  and  10 log/*0  H3H$kdk  for  re¬ 
ceiver  depths  of  450  meters  and  2500  meters,  ko  was  taken  as  5 w/c,  and 
an  effective  radius  of  100  kilometers  was  used. 
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Depth  (m)  Depth  (m) 


Wave  speed  profile  -  deep  site 


Wave  Speed  (m/s) 

Figure  6:  Profiles  of  compressions!  and  shear  wave  speeds,  Cp  and  Cs,  used  in  the  en¬ 
vironmental  model.  Compressional  wave  speeds  are  based  on  unpublished 
data  from  the  EDGE  deep  seismic  reflection  survey,  and  the  shear  wave 
speeds  were  derived  from  the  compressional  wave  speeds,  as  described  in 
the  text.  Figure  6a  is  the  deep  site  and  Fig.  6b  is  the  shallow  site. 


Bottom  Gain 


Frequency  (Hz) 


Figure  7:  Plots  of  from  data  generated  by  SAFARI  for  the  assumed  bottom 
model  at  450  m  and  2500  m  receiver  depths  for  m  =  1  and  m  =  2. 
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Figure  8:  Contours  of  the  integrands  of  Eq.  (4.41),  the  depth-dependent  Green’s 
functions  g(k,z,z'),  versus  horizontal  slowness  l/c  =  k/w  and  frequency 
for  a)  the  deep  site  (2500  m)  and  b)  the  shallow  site  (450  m). 
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Spreading  integral  versus  frequency  using  Donelan’s  empirical  model 


Figure  10:  Plot  of  10  log  Ia,t  versus  normalised  frequency  / //o  using  Donelan’s  em¬ 
pirical  model  for  the  directional  wave  spectrum  [26]  where  /0  is  double 
the  frequency  at  the  peak  in  the  wave  height  spectrum. 
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20  log  Wave  Ht.  Power  (dB)  Pred.  Noise  Power  (dB) 
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Figure  11:  Predicted  noise  level  and  inputs  for  a  2500  meter  receiver  on  1/27/91  at 
0300.  a)  Noise  spectral  level  10  logP(/)  in  dB  re  /*Pa3/Hs.  b)  Wave 
height  power  spectrum  fl(<r/2w)  in  m3/Hi.  c)  Wave  height  acoustic  input 
20logn(//2)  in  dB  re  m3/Hs.  d)  Spreading  integral  level  101og/o„(/) 
in  dB  re  Hi-1. 
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20  log  Wave  Ht  Power  (dB)  Pred.  Noise  Power  (dB) 
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Figure  12:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  1/27/91  at  0600. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  /<Pa3/H*.  b)  Wave  height  power 
spectrum  0(<y/2jr)  in  m3/Hs.  c)  Wave  height  acoustic  input  201ogR(//2) 
in  dB  re  m3/Hs.  d)  Spreading  integral  level  10 log /„„(/)  in  dB  re  Hi-1. 
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20  log  Wave  Ht.  Power  (dB)  Pred.  Noise  Power  (dB) 


Acoustic  Frequency  (Hz)  Acoustic  Frequency  (Hz) 


Figure  13:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  1/27/91  at  0900. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  /iPa3/Hs.  b)  Wave  height  power 
spectrum  f )(a/2x )  in  m3/Hs.  c)  Wave  height  acoustic  input  20  log  fl(//2) 
in  dB  re  m3/Hi.  d)  Spreading  integral  level  10  log /„„(/)  in  dB  re  Hi-1. 
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20  log  Wave  Ht.  Power  (dB)  Pred.  Noise  Power  (dB) 
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Figure  14:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  1/27/91  at  1200. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  pPa3/Hs.  b)  Wave  height  power 
spectrum  fl(<r/2ir)  in  m3/Hs.  c)  Wave  height  acoustic  input  20Iogf)(//2) 
in  dB  re  m3/Hs.  d)  Spreading  integral  level  10 log /„,,(/)  in  dB  re  Hi-1. 
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20  log  Wave  Ht.  Power  (dB)  Pred.  Noise  Power  (dB) 


Figure  15:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  3/04/91  at  0400. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  n Pa3/Hi.  b)  Wave  height  power 
spectrum  fl(<r/2ir)  in  m3/Hi.  c)  Wave  height  acoustic  input  201og0(//2) 
in  dB  re  ra3/Hi.  d)  Spreading  integral  level  10 log /„,,(/)  in  dB  re  Hi-1. 
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20  log  Wave  Ht.  Power  (dB)  Pred.  Noise  Power  (dB) 
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Figure  16:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  3/04/91  at  0600. 

a)  Noise  spectral  level  10  log  P(f)  indB  re/iPa3/Hi.  b)  Wave  height  power 
spectrum  (l(a/2x)  in  m2/Hi.  c)  Wave  height  acoustic  input  20  log  ll(//2) 
in  dB  re  iu3/Hs.  d)  Spreading  integral  level  10  log  (/)  in  dB  re  Hi-1. 
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Figure  17:  Predicted  noise  level  and  inpats  for  a  2500  m  receiver  on  3/04/91  at  0800. 

a)  Noise  spectral  level  10  log  P{f)  in  dB  re  pPa3/H*.  b)  Wave  height  power 
spectnun  fl(<r/ 2w)  in  m3/Hs.  c)  Wave  height  acoustic  input  201ogfi(//2) 
in  dB  rem3/Ht.  d)  Spreading  integral  level  10  log /„,,(/)  in  dB  re  Hi-1. 
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Figure  18:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  3/04/91  at  1000. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  fi Pa3/Hs.  b)  Wave  height  power 
spectrum  0(o/2 r)  in  m3/Hs.  c)  Wave  height  acoustic  input  201ogfi(//2) 
in  dB  re  m3/Hs.  d)  Spreading  integral  level  10  log /„„(/)  in  dB  re  H*-1. 
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Figure  19:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  3/04/91  at  1200. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  /iPa3/Hs.  b)  Wave  height  power 
spectrum  n(cr/2x)  in  m3/H*.  c)  Wave  height  acoustic  input  201ogfl(//2) 
in  dB  re  m3/Hi.  d)  Spreading  integral  level  10 log  /«„(/)  in  dB  re  Hi-1. 
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20  log  Wave  Ht.  Power  (dB)  Pred.  Noise  Power  (dB) 


Figure  20:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  3/04/91  at  1400. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  p Pa3  /Hs.  b)  Wave  height  power 
spectrum  f)(cr/2*-)  in  m3/Hs.  c)  Wave  height  acoustic  input  20  log  n(//2) 
in  dB  re  m3/Hs.  d)  Spreading  integral  level  10  log /<,„(/)  in  dB  re  Hi-1. 
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Figure  21:  Predicted  aoise  level  and  inputs  for  a  2500  m  receiver  on  3/04/91  at  1600, 
a)  Noise  spectral  level  10  log  P(f )  in  dB  re  /*Pa2/Hs.  b)  Wave  height  power 
spectrum  f 1(a/2w )  in  m2/Hi.  c)  Wave  height  acoustic  input  20 log  0(//2) 
in  dB  re  m2/Hs.  d)  Spreading  integral  level  10 log /„„(/)  in  dB  re  Hi-1. 
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Figure  22:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  2/23/91  at  0200. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  /iPa2/Hi.  b)  Wave  height  power 
spectrum  0(<7/2ir)  in  m2/Hs.  c)  Wave  height  acoustic  input  201ogn(//2) 
in  dB  re  m2/Hs.  d)  Spreading  integral  level  10 log  /<»„(/)  in  dB  re  H*-1. 
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Figure  23:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  2/23/91  at  0300. 

a)  Noise  spectral  level  10  log  P{f)  in  dB  re  ^Pa3/Hi.  b)  Wave  height  power 
spectrum  Q(<r/2*)  in  m2/Hs.  c)  Wave  height  acoustic  input  201ogfi(//2) 
in  dB  re  m2/H*.  d)  Spreading  integral  level  10  log  (/)  in  dB  re  H«-1. 
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20  log  Wave  Ht.  Power  (dB)  Pred.  Noise  Power  (dB) 
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Figure  24:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  2/23/91  at  0400. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  /iPa2/H*.  b)  Wave  height  power 
spectrum  n(<r/2jr)  in  m2/Hs.  c)  Wave  height  acoustic  input  201ogO(//2) 
in  dB  re  m2/Hs.  d)  Spreading  integral  level  10  log  /<»,,(/)  in  dB  re  Hi-1. 
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20  log  Wave  Ht.  Power  (dB)  Pred.  Noise  Power  (dB) 


Figure  25:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  2/23/91  at  0500. 

a)  Noise  spectral  level  10 log  P(/)  in  dB  re  /iPa2/Hs.  b)  Wave  height  power 
spectrum  (l(<r/2x)  in  m2/Hs.  c)  Wave  height  acoustic  input  201ogfi(//2) 
in  dB  re  m2/Hs.  d)  Spreading  integral  level  10 log  /<»„(/)  in  dB  re  Hs  *. 
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Figure  26:  Predicted  noise  level  and  inputs  for  a  2500  m  receiver  on  2/23/91  at  0600. 

a)  Noise  spectral  level  10  log  P(f)  in  dB  re  /*Pa2/Hs.  b)  Wave  height  power 
spectrum  (l(<r/2w)  in  m2/Hs.  c)  Wave  height  acoustic  input  20  log  fi(//2) 
in  dB  re  m2/Hi.  d)  Spreading  integral  level  lOlog /„,,(/)  in  dB  re  Hi-1. 
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Figure  27:  Predicted  noise  level  and  inputs  for  a  450  m  receiver  on  2/23/91  at  1000.  a) 
Noise  spectral  level  10  log  P{f)  in  dB  re  ^Pa2/H*.  b)  Wave  height  power 
spectrum  n(<r/2jr)  in  m3/Hs.  c)  Wave  height  acoustic  input  201ogfl(//2) 
in  dB  re  m2/Hs.  d)  Spreading  integral  level  10  log /<,„(/)  in  dB  re  Hi-1. 
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20  log  Wave  Ht  Power  (dB)  Pred.  Noise  Power  (dB) 


Figure  28:  Predicted  noise  level  and  inputs  for  a  450  m  receiver  on  2/23/91  at  1400.  a) 
Noise  spectral  level  10  log  P(f)  in  dB  re  jiPa2/H».  b)  Wave  height  power 
spectrum  fi(<x/2ir)  in  m2/H*.  c)  Wave  height  acoustic  input  201ogfi(//2) 
in  dB  re  m2/Hs.  d)  Spreading  integral  level  10 log  /<»„(/)  in  dB  re  Hi-1. 
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Figure  29:  Predicted  noise  level  and  inputs  for  a  450  m  receiver  on  2/27/91  at  1500.  a) 
Noise  spectral  level  10k>gP(/)  in  dB  re  /iPa2/H*.  b)  Wave  height  power 
spectrum  fl(<r/2ir)  in  m2/Hi.  c)  Wave  height  acoustic  input  20  log  0(//2) 
in  dB  re  m2/Hi.  d)  Spreading  integral  level  10 log  /<,„  (/)  in  dB  re  Hs-1. 
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Figure  30:  Predicted  noise  level  and  inputs  for  a  450  m  receiver  on  2/27/91  at  1800.  a) 
Noise  spectral  level  101ogP(/)  in  dB  re  ^Pa3/H«.  b)  Wave  height  power 
spectrum  fl(<r/2ir)  in  m3/Hx.  c)  Wave  height  acoustic  input  201ogn(//2) 
in  dB  re  m3/Hs.  d)  Spreading  integral  level  10 log  /„„  (/)  in  dB  re  H*_1. 
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Wind  speed  and  wave  height,  1/27/91,  Discus  E 


Figure  31:  Meteorological  data  as  measured  at  Discus  E  on  1/27/91.  a)  Wind  vector. 

The  direction  of  the  arrow  indicates  the  direction  of  the  wind,  with  a 
downward  arrow  indicating  a  wind  from  the  west.  The  magnitude  of  the 
arrow  indicates  the  wind  speed,  with  the  wind  speed  scale  given  on  the 
vertical  axis,  b)  Wind  speed  and  wave  height. 
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Figure  32:  Meteorological  data  as  measured  at  Discus  E  on  3/04/91.  a)  Wind  vector. 

The  direction  of  the  arrow  indicates  the  direction  of  the  wind,  with  a 
downward  arrow  indicating  a  wind  from  the  west.  The  magnitude  of  the 
arrow  indicates  the  wind  speed,  with  the  wind  speed  scale  given  on  the 
vertical  axis,  b)  Wind  speed  and  wave  height. 


92 


Significant  Wave  Height  (m) 


Wind  Speed  (m/s)  Wind  Speed  (m/s) 


Wind  vector,  1200  2/22/91  - 1200  2/23/91,  Discus  E 


Time  of  Day 


Wind  speed  and  wave  height,  1200  2/22/91  -  1200  2/23/91,  Discus  E 


Tune  of  Day 

Figure  33:  Meteorological  data  as  measured  at  Discus  E  on  2/22/91*2/23/91.  a) 
Wind  vector.  The  direction  of  the  arrow  indicates  the  direction  of  the 
wind,  with  a  downward  arrow  indicating  a  wind  from  the  west.  The 
magnitude  of  the  arrow  indicates  the  wind  speed,  with  the  wind  speed 
scale  given  on  the  vertical  axis,  b)  Wind  speed  and  wave  height. 
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Significant  Wave  Height  (m) 


ONR  OBS  combined  filter  response 


Figure  34:  Amplitude  response  of  the  anti-aliasing  and  pre-whitening  filter  used  in 
the  ECONOMEX  instruments  with  an  8  Hs  sampling  rate.  Response 
amplitude  in  dB  re  1  Volt. 
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Noise  Spectral  Density  (dB) 


Estimated  spectrum  OBS  58  DPG  3/4/91  1700 


Figure  35:  An  example  of  a  full  spectral  estimate  using  data  from  OBS  58  DPG  on 
3/04/91  at  about  1700  hours.  This  spectrum  was  generated  as  described 
in  the  text  with  the  exception  that  a  2048  point  FFT  was  used  versus  a 
512  point  FFT.  The  spectral  level  is  in  dB  re  pPa2/Hs. 
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Noise  Power  (dB) 


2500  m,  1/27/91, 1800 


450  m,  2/23/91, 1400 


Frequency  (Hz) 

Frequency  (Hz) 

Figure  36:  Predicted  and  observed  noise  lev¬ 
els,  1/27/91,  1800,  OBS  58  (2500 
m),  in  dB  re  #iPa2/H*.  Observed 
spectrum  is  solid  curve. 

Figure  37:  Predicted  and  observed  noise  lev¬ 
els,  2/23/91,  1400,  OBS  63  (450 
m),  in  dB  re  ^Pa2/H*.  Observed 
spectrum  is  solid  curve. 

Noise  Power  (dB) 


2500  m,  2/23/91,0700 


Frequency  (Hz) 


Figure  38:  Predicted  and  observed  noise  lev¬ 
els,  2/23/91,  0700,  OBS  58  (2500 
m),  in  dB  re  jiPa2/Hs.  Observed 
spectrum  is  solid  curve. 


2500  m,  1/27/91, 1200 


Figure  39:  Predicted  and  observed  noise  lev¬ 
els,  1/27/91,  1200,  OBS  58  (2500 
m),  in  dB  re  /<Pa2/Hi.  Observed 
spectrum  is  solid  curve. 
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Noise  Power  (dB) 


2500  m,  3/04/91, 0000 


Figure  40:  Predicted  and  observed  noise  lev¬ 
els,  3/04/91,  0000,  OBS  58  (2500 
m),  in  dB  re  /xPa2/H*.  Observed 
spectrum  is  solid  curve. 


450  m,  2/27/91, 1800 


Frequency  (Hz) 


Figure  41:  Predicted  and  observed  noise  lev¬ 
els,  2/27/91,  1800,  OBS  63  (450 
m),  in  dB  re  /*Pa2/Hi.  Observed 
spectrum  is  solid  curve. 
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Noise  Power  (dB) 


Noise  Power  (dB) 


2500  m,  1/27/91,0300 


Figure  44:  Predicted  and  observed  noise  lev¬ 
els,  1/27/91,  0300,  OBS  58  (2500 
m),  in  dB  re  jrPa2/Hi.  Observed 
spectrum  is  solid  curve. 


2500  m,  1/27/91,0600 


Figure  45:  Predicted  and  observed  noise  lev¬ 
els,  1/27/91,  0600,  OBS  58  (2500 
m),  in  dB  re  /rPa2/F*.  O^servAd 
spectrum  is  solid  curve. 
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2500  m,  1/27/91,0900 


Frequency  (Hz) 


Figure  46:  Predicted  and  observed  aoise  lev¬ 
els,  1/27/91,  0900,  OBS  58  (2500 
m),  in  dB  re  «Ps3/Hr.  Ofcccrvcd 
spectrum  is  solid  curve. 


2500  m,  3/04/91, 0600 


Frequency  (Hz) 


Figure  48:  Predicted  and  observed  noise  lev¬ 
els,  3/04/91,  0600,  OBS  58  (2500 
m),  in  dB  re  /*Pa2/Hi.  Observed 
spectrum  ia  solid  curve. 


2500  m,  3/04/91, 0800 


Figure  49:  Predicted  and  observed  noise  lev¬ 
els,  3/04/91,  0800,  OBS  58  (2500 
m),  in  dB  re  j»Pa2/Hi.  Observed 
spectrum  is  solid  curve. 


102 


Noise  Power  (dB) 


2500  m,  3/04/91, 1000 


Figure  50:  Predicted  and  observed  noise  lev¬ 
els,  3/04/91,  1000,  OBS  58  (2500 
m),  in  dB  re  /*Pa2/Hs.  Observed 
spectrum  is  solid  curve. 


2500  m,  2/23/91, 0200 


Frequency  (Hz) 


Figure  51:  Predicted  and  observed  noise  lev¬ 
els,  2/23/91,  0200,  OBS  58  (2500 
m),  in  dB  re  /«Pa2/Hi.  Observed 
spectrum  is  solid  curve. 


Noise  Power  (dB) 


2500  m,  2/23/91,0300 


Frequency  (Hz) 


Figure  52:  Predicted  and  observed  noise  lev¬ 
els,  2/23/91,  0300,  OBS  58  (2500 
m),  in  dB  re  #xPa2/Hs.  Observed 
spectrum  is  solid  curve. 


2500  m,  2/23/91,0400 


Figure  53:  Predicted  and  observed  noise  lev¬ 
els,  2/23/91,  0400,  OBS  58  (2500 
m),  in  dB  re  /<Pa2/Hs.  Observed 
spectrum  is  solid  curve. 
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Noise  Power  (dB) 


Noise  Power  (dB) 


